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GENERALIDADES 


Desde los primeros pasos en el cálculo diferencial, de todos es conocido que, dada una 
función y — f(x), su derivada = f'(x) es también una función que se puede encontrar 

mediante ciertas reglas. Por ejemplo, si y = e~ x , entonces = —3x 2 e~ x o, lo que es 
lo mismo, 4 ^ = —3 x 2 y. El problema al que nos enfrentamos ahora no es el de calcular 
derivadas de funciones; más bien, el problema consiste en: si se da una ecuación como 
= —3 x 2 y, hallar de alguna manera una función y = f(x) que satisfaga dicha ecuación. 
En una palabra, se desea resolver ecuaciones diferenciales. 

La forma de ecuación diferencial más sencilla que puede pensarse es ^ = f(x). 
Resolverla consiste en encontrar una función cuya derivada sea f(x). es decir, encontrar 
las primitivas (integrales indefinidas) de f(x). Por tanto, podemos decir que los métodos 
de resolución de ecuaciones diferenciales constituyen una generalización del cálculo de 
primitivas. 

Definición 1. Llamamos ecuación diferencial (E. D.) a una ecuación que relaciona 
una función (o variable dependiente), su variable o variables (variables independientes), y 
sus derivadas. Si la ecuación contiene derivadas respecto a una sola variable independiente 
entonces se dice que es una ecuación diferencial ordinaria (E. D. O.); y si contiene las 
derivadas parciales respecto a dos o más variables independientes se llama ecuación en 
derivadas parciales (E. D. P.). 


Ejemplos de ecuaciones diferenciales ordinarias son 

^ — 4y = 2, (x + 2 y) dx — 3y dy = 0 

dx 


y 


mientras que 


y 


fy 

dx 2 


dy 

dx 


+ 3y = 0; 


du du 

x-r r + y^r = u 
dx dy 


d 3 u d 2 u du 
dx 3 dt 2 dt 


son ecuaciones en derivadas parciales. 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 
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Métodos clásicos de resolución de E. D. O. 


Otro tipo de ecuaciones que pueden estudiarse son las ecuaciones diferenciales de 
retraso (o retardo ), como es el caso de 

u'(t) = 7-2u(t-3). 

Están caracterizadas por la presencia de un desplazamiento t — tg en el argumento de la 
función incógnita u(t). En general, son más difíciles de manejar que las E. D. sin retraso. 
No nos ocuparemos aquí de ellas. 

Definición 2. Se llama orden de la ecuación diferencial al orden de la derivada o derivada 
parcial más alta que aparece en la ecuación. 

Así, por ejemplo, las ecuaciones (1) y (3) son de orden 1, (2) es de orden 2 y (4) de 
orden 3. 

En lo que sigue nos preocuparemos sólo de ecuaciones diferenciales ordinarias y, como 
no habrá lugar a confusión, las denominaremos simplemente E. D. Por lo general, salvo 
que el contexto nos indique otra notación (o ésta provenga de los cambios de variable que 
efectuemos), utilizaremos x para denotar la variable independiente e y para la variable 
dependiente. 

Definición 3. Decimos que una ecuación diferencial (de orden n ) está expresada en forma 
implícita cuando tiene la forma 


F(x,y,y',...,y^) = 0 

siendo F una función F:Q C R n+2 — y M con O un subconjunto (generalmente abierto) 
de M n+2 . Y decimos que está expresada en forma explícita cuando tenemos 

y in) = f(x,y,y',..., j/ ( “- 1) ) 


con f:D C R n+1 — y M una función deñnida en un subconjunto D (generalmente abierto) 
de R n+1 . 


Una clase importante de E. D., bien estudiada y con buenas propiedades, es la 
siguiente: 


Definición 4. Se dice que una ecuación diferencial es lineal si tiene la forma 


, . d n y . . d n 1 y 

an{ - X) dx¿ + an - l{ - X) dD^ I + 


+ a\(x)^j- + a 0 (x)y = g(x); 
dx 


y se llama lineal homogénea si, además, g(x) = 0. Dada una ecuación lineal, su 
correspondiente ecuación lineal homogénea en la que se ha hecho g(x) — 0 se denomina 
lineal homogénea asociada. Una ecuación que no es lineal se dice no lineal. 


Nuestro objetivo es resolver ecuaciones diferenciales, esto es, encontrar sus soluciones. 




Generalidades 
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Definición 5. Decimos que una función y = ip(x) definida en un intervalo I (es decir, 
ip: I C R — > R ) es solución de una ecuación diferencial en el intervalo si, sustituida en 
dicha ecuación, la reduce a una identidad. (En otras palabras, si satisface la E. D.) Una 
E. D. se dice resoluble (o integrable) por cuadraturas si su solución es expresable 
mediante integrales. 

En general, la solución de una ecuación diferencial de orden n dependerá de n 
parámetros. Pero incluso de esta forma pueden no obtenerse todas las soluciones de una 
E. D. Por ejemplo, cuando tenemos una familia uniparamétrica de soluciones de una E. D., 
una sencilla interpretación geométrica nos muestra que también la envolvente de la familia 
de curvas (si existe) es solución de la E. D. 

A continuación, nos dedicaremos a explicar los diversos métodos clásicos de resolución 
de E. D. No liaremos hincapié en el intervalo de definición de las soluciones, ni efectuaremos 
un estudio detallado de la rigurosidad de los métodos empleados que, en esencia, descansan 
siempre en la regla de la cadena y los teoremas de la función inversa y de la función 
implícita. No nos detendremos nunca en comprobar las hipótesis de estos teoremas, sino 
que supondremos en todo momento que las funciones que aparecen en los métodos descritos 
son lo suficientemente «buenas», o están lo suficientemente restringidas en su dominio, para 
que siempre se satisfagan las hipótesis necesarias. Tampoco nos preocuparemos en exceso 
de saber si hemos obtenido todas las soluciones; a este respecto, en algunos casos nos 
interesaremos por las soluciones singulares de una E. D., como puede ser la envolvente de 
una familia de soluciones. 

Nos apresuramos a señalar que las fórmulas generales que aparecen como solución 
de diversos tipos de ecuaciones no deben memorizarse; más bien, el procedimiento debe 
desarrollarse completo cada vez. Para ello bastará recordar unos cuantos puntos esenciales 
que destacamos en las cajas de texto que hemos denominado «recetas». 

Por último, comentar que, en los ejemplos que nos aparecerán, el lector puede 
entretenerse en representar gráficamente las soluciones de las E. D. planteadas, al menos 
en los casos más sencillos o efectuando un simple bosquejo de su apariencia. No pensemos 
en esto como una pérdida de tiempo, pues ayuda a comprender la naturaleza del problema 
y de sus soluciones. 





ECUACIONES EXPLÍCITAS DE 
PRIMER ORDEN 


Son las que tienen la forma 


y' = f(x,y). 


APARTADO 1. 

Variables separadas. 

Si tenemos la E. D. 

g{x) = h(y)y ', 

formalmente, podemos poner g(x) dx = h(y) dy\ si suponemos que G es una primitiva de 
g y H una de h, tendremos G'(x ) dx — H'(y ) dy e, integrando, G(x) — H(y) + (7, que es 
la solución general de la ecuación. 

Expliquemos con un poco más de rigor por qué funciona el método: Sea y — p(x) 
una solución de la E. D., es decir, p(x) debe cumplir g(x) — h((p(x))<p'(x). Pero H es 
una primitiva de h. así que, por la regla de la cadena, g(x) — h(p(x))p' (x) = ( H o tp)'(x). 
Integrando, G(x) = ( H o ip){x) + C (lo que antes hemos expresado como G(x) — H(y)+C), 
de donde <p{x) = H~ 1 (G(x) — C ). 

En los pasos anteriores, está justificado emplear la regla de la cadena cuando ip y H 
son derivables, lo cual es cierto sin más que suponer que h sea continua. Y finalmente, para 
poder despejar ip mediante el uso de II ~ 1 bastaría con exigir además que h no se anulara 
en el intervalo de definición, con lo cual, como H r = h ^ 0, H es creciente o decreciente 
luego existe H -1 (en otras palabras, como la derivada de H no se anula, el teorema de la 
función inversa nos asegura que existe II 1 )- 

Las ecuaciones en variables separadas son las más sencillas de integrar y, a la vez, 
las más importantes, ya que cualquier otro método de resolución se basa esencialmente en 
aplicar diversos trucos para llegar a una ecuación en variables separadas. En ellas hemos 
visto, con todo rigor, qué hipótesis hay que imponer para que el método que conduce 
a la solución esté correctamente empleado, y cómo se justifica el funcionamiento del 
proceso. A partir de ahora no incidiremos más en estos detalles que, aunque importantes, 
sobrecargarían la explicación. El lector puede detenerse mentalmente a pensar en ellos, 
justificando adecuadamente los pasos que se efectúen. 


5 



6 


Métodos clásicos de resolución de E. D. O. 


En cualquier caso, conviene recordar que la expresión 4^ es simplemente una útil 
notación para designar la derivada de y respecto de x, no un cociente de dy dividido 
por dx : ni dy ni dx tienen entidad en sí mismas. Esta notación se emplea, no porque 
sí ni para introducir confusión, sino que, al contrario, se usa porque es consecuente con 
los enunciados de varios importantes resultados. Ya Iremos visto cómo resulta adecuada 
a la hora de recordar cómo resolver ecuaciones en variables separadas g(x) — h(y') 
descomponiendo g(x) dx — h(x) dy (como si fuese realmente una fracción) e integrando 
ambos lados de la expresión anterior. Pero no sólo aquí se pone de manifiesto la utilidad 
de esta notación. Por ejemplo, el teorema de la función inversa prueba (con las hipótesis 
adecuadas) que, cuando y es una función de x, si se despeja x como función de y se cumple 

dx 1 1 

dy dy / dx y'{x)' 

es decir, se produce un comportamiento similar a si estuviéramos operando con fracciones. 
Análogamente, si tenemos que z es una función de y y, a su vez, y una función de x. la 
regla de la cadena establece que la derivada de la función compuesta z(x) es 

dz dz dy 
dx dy dx ’ 

que es como si simplificáramos dy en los supuestos cocientes de la derecha. Esto permite 
usar las notaciones del tipo -p y su comportamiento como si fuesen fracciones como regla 
nemotécnica de los resultados anteriores. 



RECETA 1. Variables separadas. 

Son de la forma 

g{x) = h{y)y' . 

Formalmente, se separa g(x ) = h(y) -A en g{pc) dx — h{y) dy y se 
integra. 


Ejemplo 1. Resolver 

— + (sens)y = 0. 
dx 

Despejando, ^ = —(sena;) dx e, integrando, log y — cosa; + C, es decir, y = e cosx + G . 
Sin más que tomar K = e c encontramos las soluciones y = Ke cosx . Fijarse que, en 
principio, parece que K tiene que ser positiva; pero en realidad la integral de — es log \y\, 
lo que nos llevaría a soluciones con valores negativos de K. Por último, notar y = 0 (es 
decir, tomar K = 0) también es claramente una solución de la E. D., aunque no se obtiene 
con el método seguido. Así pues, la solución general de la E. D. es de la forma y — Ke cosx 
con K G IR. 





E. D. explícitas de primer orden 
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APARTADO 2. 

Ecuación de la forma y' — f(ax + by). 

Si a = 0 o b — 0, la ecuación es separable. En otro caso, efectuemos el cambio de 
función y(x) por z(x) dado por z — ax + by , de donde z' — a + by' y, por tanto, y ' = a 7 a . 
Entonces, sustituyendo en la E. D. obtenemos = f(z), es decir, z' — a + bf(z ), que 
es de variables separadas. La escribimos como 


dz 

a + bf(z) J 

con lo que, integrando, x — f (a + bf(z))~ 1 dz — (¡>{z,C). Así pues, las soluciones de la 
E. D. de partida serán 

x — (f>(ax + by, C ), 

de modo que hemos encontrado y como función de x expresada en forma implícita. 


RECETA 2. Ecuación de la forma y 1 = f(ax + by). 

El cambio de función y(x) por z(x) dado por z — ax + by la 
transforma en una de variables separadas. 


Ejemplo 2. Resolver 

y' -e x e y = -1. 

Tenemos y' + 1 = e x+y , con lo que si efectuamos el cambio de función dado por la 
sustitución z = x + y, la ecuación queda transformada en z' — e z , es decir, dx = e~ z dz, 
ecuación en variables separadas cuya solución es x — —e~ z + C. Volviendo a las variables 
iniciales, C — x = e~ x ~ y , de donde log(C — x) = —x — y, y por tanto la solución de la 
E. D. de partida es y — — log(C — x) — x. (Observar que no nos hemos preocupado — ni lo 
haremos de aquí en adelante — de poner módulos cuando al calcular una integral aparece 
un logaritmo. El lector podría analizar estos casos con mucho más cuidado.) 


APARTADO 3. 

Homogéneas. 

Supongamos que tenemos la ecuación 





Métodos clásicos de resolución de E. D. O. 


Para resolverla, hacemos el cambio de función y(x) por u(x ) mediante u = Así, derivando 
y — ux tenemos y' = u'x + u, es decir, u'x + u = f(u). Esta ecuación, que podemos poner 
como u'x — f(u) — u, es de variables separadas. Vamos a solucionarla: 

• Si f(u) ^ u, podemos escribir f $_ u = f e, integrando, f f ¿)- u = lo g(#)- 
Despejando x obtenemos x = Ce^ u ^ con <j)(u) = f f^)- u • P° r t an t°: las curvas con 
ecuaciones paramétricas 

j x = 

\y = Cue^ u) 

son solución de la ecuación diferencial para cada C G I. (Esto constituye una familia 
de curvas liomotéticas: una curva se obtiene de otra mediante una liomotecia, es decir, 
multiplicando los valores de x e y por una constante.) A veces, es conveniente expresar 
estas soluciones de otras formas. Siempre puede ponerse x = Ce^ y / X \ solución dada 
mediante una función implícita. Y, cuando en x — Ce+M se logra despejar de alguna 
forma u — H(x,C), la solución de la E. D. queda mucho más sencilla: y — xH(x, C ). 

• Supongamos ahora que existe algún u$ tal que f(uo) — u$. En este caso, es inmediato 
comprobar que la recta y = uqx es solución: y' — uq = /(«o) = /(f ), luego se satisface 
la ecuación diferencial. Este tipo de soluciones que no se obtienen con el procedimiento 
general suelen denominarse soluciones singulares. 

Nota: En general, una función h(x,y ) se dice homogénea de grado a si h(Xx,Xy) — 
X a h(x,y). Es inmediato comprobar que una E. D. de la forma 

P(x, y) dx + Q{x, y) dy — 0 

con P(x, y) y Q(x, y) funciones homogéneas del mismo grado es, efectivamente, una 
ecuación diferencial homogénea (despejar y' — ^ — — q^’x^/x) ) Y extraer 

A = x de P y Q). De aquí proviene el nombre de este tipo de ecuaciones. 


RECETA 3. Homogéneas. 

Son de la forma 



Se hace el cambio de función y(x) por u(x) mediante y 
transformándose así la E. D. en una de variables separadas. 


ux, 


/ 

y 


2 xy-y 


ar 


2 


Ejemplo 3. Resolver 
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Con el cambio y — ux podemos poner y' = 2^ — (^) 2 = 2u — u 2 . Como y' — u'x + u, 
sustituyendo tenemos u'x + u — 2u — u 2 , es decir, xu' — u — u 2 . 

Si u ^ u 2 , podemos poner Para integrar, descomponemos j 1 2 = ^ + 

lo que se satisface para A — B = 1. Entonces, integrando, log u — log(l — u) = log es 

decir, y sustituyendo u — ^ tenemos i^!y/ x — fj, de donde Cy = x{x — y). De 

aquí es fácil despejar explícitamente y si así se desea. 

Por otra parte, a partir de Uq = 0 y Uq — 1 (para las cuales u = u 2 ), se tienen las 
soluciones singulares y = 0 e y — x. 


APARTADO 3'. 

Reducibles a homogéneas. 

Consideremos la ecuación 


, _ . í g\x + hy + Ci \ 

\ ax + by + c ) 

Para resolverla, hay que distinguir dos casos: 

3'.1. Supongamos en primer lugar que las rectas ax + by + c — 0 y a±x + b\y + C\ — 0 
se cortan en el punto (xo, yo)- Así, tendremos que ax + by + c = a(x — xo) + b(y — yo) y 
a\x + b\y + c\ — a\{x — xo) + b\ (y — yo)- Hagamos ahora el cambio de variable y de función 
X = x — Xq , Y — y — yo, con lo cual 


y, _ / _ r ( ai(x - Xq) +bi(y - y 0 ) \ = í aiX + btY \ / + 

V a ( x - x o) + Hy - yo) ) \ aX + bY ) y a + bY J ’ 

es decir, hemos reducido la ecuación a una homogénea. 

3'. 2. En segundo lugar, supongamos que ax + by + c — 0 y aiX + biy + Ci = 0 son rectas 
paralelas, con lo cual podrá ponerse (ai, b\ ) = K(a, b) para algún K e M. Efectuemos ahora 
el cambio de función z — ax + by. Derivando, z' = a + by' , o sea, y' = 2 7 a . Si sustituimos 
en la E. D. original obtenemos 


dz 

dx 


a + bf 


Kz + ci 


que es de variables separadas. 


z + c 
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Métodos clásicos de resolución de E. D. O. 


RECETA 3'. Reducibles a homogéneas. 

Son de la forma 

/ _ , f aix + hy + ci \ 

\ ax + by + c ) 

3'.1. Si las rectas ax + by + c — 0 y a\x + b\y + c\ — 0 se cortan 
en (xo,yo), se hace el cambio de variable y de función X = x — xo, 
Y = y — y 0 . La ecuación se reduce a una homogénea. 

3'. 2. Si ax + by + c — 0 y aix + biy + ci — 0 son rectas paralelas, se 
hace el cambio de función z — ax + by. La nueva ecuación que aparece 
es de variables separadas. 


Ejemplo 3'.1. Resolver 


— 2x + Ay — 6 
x + y — 3 


Las rectas —2x + Ay — 6 — 0 y x + y — 3 = 0 se cortan en el punto (x,y) — (1,2), con 
lo que efectuamos el cambio X = x — 1, Y = y — 2, Y' = y' . Sustituyendo, obtenemos la 
ecuación homogénea 

, -2X + 4 Y -2 + AY/X 

~ X + Y ~ 1 + Y/X ' 

Para resolverla, hacemos un nuevo cambio u — Y / X . de donde Y = aX . Y' = u'X + u. 
Tras sustituir, tenemos u'X + u — que, a la postre, podemos poner como —X ^ = 

m 2 -3m+2 
u + 1 

Tenemos ahora que distinguir cuándo, y cuándo no, se anula la expresión u 2 — 3u + 2. 
Resolviendo u 2 — 3u + 2 = 0, esto ocurre para u — 1 y u — 2. 

Analicemos en primer lugar el caso u 2 — 3u + 2 ^ 0. Así, podemos escribir 


dX 

~X 


— ( u + 1) du 
u 2 — 3u + 2 


A 


u- 1 


du + 


B 

u — 2 


du = 


2 


u - 1 


du — 


3 

u — 2 


du 


de donde, integrando, 21og( , u — 1) — 31og(w — 2) = log(JiA) y, consiguientemente, 
| ^ ~ 2 j 3 = KX. Sustituyendo ahora u = Y/X llegamos fácilmente a (Y— X) 2 = K(Y— 2X) 3 ; 
y volviendo a las variables originales x e y obtenemos las soluciones (y—x— l) 2 = K(y— 2x) 3 
de la E. D. de partida. 

Finalmente, con uq — 1 y uq — 2 tenemos, respectivamente, las soluciones Y = X e 
y = 2X que, sustituyendo X e Y por su valor, se traducen en y = x + 1 y y — 2x. 
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Ejemplo 3'. 2. Resolver 

/ x-y - 1 

V = o- 

x — y — 2 

Efectuamos el cambio de función z — x — y, de donde z' = 1 — y' . Sustituyendo, 
tenemos —z' + 1 = y consiguientemente — — 1 = , que es la ecuación 

(z — 2) dz — —dx, cuyas variables están separadas. Integrando, ^(z — 2) 2 = —x + K . y 
finalmente, sustituyendo de nuevo z — x — y y denotando C — 2 K, obtenemos que las 
soluciones de la E. D. original son (x — y — 2) 2 + 2x — C. 


APARTADO 3". 


Homogéneas implícitas. 

Sea la ecuación 



Para resolverla, consideremos la curva F(ot,¡3) = 0 y supongamos que hemos logrado 
encontrar una representación paramétrica de la curva dada por a = <p(t), (3 = ip{t). Es 
decir, que se verifica F(<p(t),ip(t)) = 0. Hagamos ahora el cambio de función y por t 
mediante y = <p(t), teniendo en cuenta que y 1 — ^{t). 

Si derivamos y — xip(t) respecto de x tenemos y' — < p(t ) + x(p\t)4^, es decir, 
V’(í) = <p(t) + x<p , (t)^, o 

(jf 

^(í) - <p(t) = X(p'(t] ) — 


que, en principio, es una ecuación en variables separadas. 

• Si ip(t) ^ <p(f), podemos poner = yy, cuya solución será x = con 

(f)(t) = f • De aquí que la E. D. de partida tiene las soluciones 


x = Ce*® 

y = C(p{t)e m ’ 


• Si existe t 0 tal que ip(to) = <p(fo), formalmente, podemos pensar 0 = x(p'(t)-^, luego 
<p'(t) — 0 y por tanto <p(t) — cte. = <p(to), lo que nos llevaría a la solución y = x<p(to). Esta 
recta es, efectivamente, una solución de la E. D., como podemos comprobar directamente: 
F (x’V') = F (<P(to),<p(to)) = F(<p(t 0 ),i/>(t 0 )) = 0 . 
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Métodos clásicos de resolución de E. D. O. 


RECETA 3". Homogéneas implícitas. 

Son de la forma 



Consideramos la curva F(oe, d) = 0. Si encontramos una representación 
paramétrica a = <p(t), ¡3 — F(<p(t),ip(t)) = 0, se hace el cambio de 

función y por t mediante ^ = <p(t), y' — ^(t). Así, derivando y — x(p(t) 
respecto de x, aparece una ecuación en variables separadas. 


Ejemplo 3" . Resolver 

x 2 (y') 2 -(y 2 + x2 ) = o. 

Si ponemos la ecuación en la forma ( y ’) 2 — (^) 2 = 1, podemos recordar que el coseno 
y el seno hiperbólicos satisfacen la relación (chí) 2 — (shf) 2 = 1, lo que se adecúa a nuestras 
necesidades. Así, tomemos ahora - = shí y y' — cli t. Si derivamos y — xsht tenemos 
y' — shí + xdht 4^, o sea, chf = shí + xcht Despejando, ^ = ch t dt (notar que, 
en esta ecuación, el denominador no se anula nunca ya que chí > sil t. luego no hay que 
preocuparse de analizar por separado las raíces de ch t — sh t, — 0) e, integrando, x = Ce r - L ) 
con 

m= l dTr=lh7* = / 1 ¿5-* = i/(l + e 2 ‘)* = it+ie“ 

Entonces, la E. D. original tiene como soluciones las curvas 

íx = Ce i / 2+e2t / 4 
\ y = C(sht)e t/2+e2i/4 ' 


APARTADO 3"' . 

Sea la ecuación y 1 = f(x, y) 

con / tal que, para algún a fijo, verifica 

f(\x,\ a y) = X a_ 1 /(a;,y). 

Nótese en primer lugar que, cuando a — 0, sin más que tomar X — x tenemos y' = f(x. y) = 
x~ x f(l,y), que es una ecuación en variables separadas; y, cuando o: = 1, se puede poner 

y' = f(x, y) = f (x, x^j = f (l, , 

es decir, nos encontramos ante una E. D. homogénea. En otro caso, veamos cómo el cambio 
de función y — z a transforma la ecuación en una homogénea: 
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Derivando, y' = az a l z' y, sustituyendo en la E. D. original, az a x z' — f(x,z a ), es 
decir, 

~ (~) a l f( x i z<x ) = -f (-x, (-) V) = -f (-, l) 

ql\zJ a \z \zJ J a \z / 

que, efectivamente, es una E. D. homogénea. 

Lógicamente, como, al hacer en la homogénea el cambio z — ux, ésta se transforma en 
una de variables separadas, si hubiéramos efectuado desde el principio el cambio y = (ux) a , 
nuestra E. D. se hubiera convertido directamente en una de variables separadas. 

Por último comentar que, extrayendo A = x, este tipo de ecuaciones puede ponerse 
como 

y' = f (x, y) = f (x, x a ^ = X a - 1 f (l, - x a ~ l h ( . 

Pero, al menos a simple vista, no parece más sencillo describir las ecuaciones que estamos 
tratando como «las que tienen la forma y ' = x a ~ 1 h(yx~ a )» en lugar de como lo liemos 
hecho en este apartado. 


RECETA 3"'. Si la ecuación y' — f(x,y) 
es tal que, para algún a ^ 0 fijo, / satisface 

f(Xx,X a y) = 

entonces el cambio de función y — z a transforma la ecuación en una 
homogénea. (Si a = 1, la E. D. ya es homogénea; y si / cumple la 
relación anterior con a — 0, la E. D. es de variables separadas.) 


Ejemplo 3'". Resolver 


/ 

y 


1 y _ 0 Vx 

2 x y 2 


Dada f(x. y) — \x l y — ‘ix 1 l 2 y 2 , para intentar encontrar a tanteamos /( Xx, X a y) = 
7¿(Xx)~ 1 (X a y) — 3(Ax) 1 / 2 (A a y)~ 2 = ^X a ~ 1 x~ 1 y — 3X 1 ^ 2 ^ 2a x 1 ^ 2 y~ 2 , y observamos que esto 
es igual a A “ _1 /(x, y) sin más que tomar a = 

Entonces, si sustituimos y — z 1 ^ 2 en la E. D., tenemos \z~ 1 ! 2 z' = \x~ x z x ! 2 — 
3x x I 2 z ~ x , es decir, z' — ^ — 6(^) -1 / 2 , que es homogénea. Para resolverla, tomamos ahora 
z = ux, con lo cual, sustituyendo, u'x+u = u— 6n -1 / 2 , o sea, u 1 ^ 2 du — — 6x _1 dx, ecuación 
en variables separadas. Integrándola, | u 3//2 — — 61og(^), luego x = Cexp(— w 3 / 2 /9). 
Deshaciendo los cambios z = ux y y — ux encontramos que las soluciones de la E. D. 
de partida son, expresadas como curvas en paramétricas, 

i x = Ce~ u3/2 / 9 

\ y = ' 
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Métodos clásicos de resolución de E. D. O. 


APARTADO 4. 

Ecuaciones exactas. 

Llamamos exacta a una ecuación diferencial 

P(x, y) dx + Q(x, y) dy = 0, 

es decir, y' = ^ , que cumple Py = Q x (con la notación P y = ^,Q X = |g). 

Antes de explicar cómo resolverlas, comentemos brevemente algo sobre «expresiones 
diferenciales» (rigurosamente hablando, estamos tratando con 1 -formas diferenciales w = 
P dx + Q dy, aunque no entraremos en ello): 

Supongamos de antemano en todo lo que sigue que P y Q son de clase C 1 (continuas 
con derivadas parciales continuas) en su dominio de definición (un abierto de M 2 ). Una 
expresión diferencial P(x, y) dx + Q(x, y) dy se dice que es una diferencial cerrada en una 
región R del plano xy si se verifica P y [x,y) — Q x (x,y ) para todo (x,y) & R. Y se dice 
exacta en R cuando existe alguna función F(x,y ) tal que ^ = P y QP = Q para todo 
(x,y) G R; en otras palabras, si la diferencial de F es dF = P dx + Q dy (. F , que es 
única salvo constantes, se denomina función potencial). El teorema de Schwartz sobre 
igualdad de derivadas cruzadas nos asegura que cualquier expresión diferencial exacta 
es cerrada. Lo contrario no es cierto en general, aunque sí en una clase muy amplia 
de dominios de IR 2 : los simplemente conexos que, intuitivamente, son los que no tienen 
agujeros. Demostrar este hecho no es excesivamente sencillo, pero tampoco es necesario 
para lo que aquí pretendemos. En realidad, el lema de Poincaré (que normalmente se 
prueba en cualquier curso de cálculo integral en varias variables) asegura que una expresión 
cerrada es exacta siempre que el dominio sea estrellado , lo que significa que exista un 
punto del dominio que se pueda unir a todos los demás mediante un segmento sin salimos 
del dominio; en particular, los conjuntos convexos son estrellados. Además, esto asegura 
que, dada cualquier expresión cerrada, es exacta localmente, es decir, alrededor de cada 
punto podemos restringir el dominio de tal forma que la expresión sea exacta en ese nuevo 
dominio. Por lo tanto, en lo que a nosotros concierne, podemos identificar los conceptos 
de exacto y cerrado, ya que no nos estamos preocupando de dónde están definidas las 

E. D. que tratamos de resolver ni en qué intervalo existen las soluciones. En realidad, en 
ecuaciones diferenciales suele hablarse siempre de exacto aún refiriéndose a que se satisface 
la igualdad P y = Q x . 

Una E. D. exacta es una expresión exacta igualada a cero. Veamos cómo resolverlas: 
si tenemos P dx + Q dy = 0 exacta, como existe F tal que dF = P dx + Q dy, entonces la 
ecuación podemos ponerla en la forma dF = 0 y, por tanto, su solución será F(x,y) — C 
(siendo C constante arbitraria). Así pues, basta con que encontremos la función potencial 

F. El procedimiento para hallarla que, según veremos, funciona gracias a que P y — Q x , es 
como sigue: 

Buscamos F tal que QF = P; así, es posible encontrar F integrando P(x,y) respecto 
a x mientras se mantiene y constante, es decir, F(x,y) = f P(x,y) dx + <p(y), donde la 
función arbitraria <p(y) es la «constante» de integración. Derivando respecto de y obtenemos 
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= ~§y f ^( X ’y) d x + p'iy)- P° r °t ra parte, si utilizamos ^ = Q, de aquí resulta que 
<p'(y) — Q{x,y) — f P(x,y ) dx; ésta es realmente una expresión independiente de x ya 
que 


d_ 

dx 




dQ d / d 
dx dy V dx 



dQ__dP_ 
dx dy 


Una vez conocida <p'(y), integrando obtenemos <p(y) y, sustituyendo su valor, llegamos a 
la función potencial F(x,y). Así, quedan halladas completamente las soluciones buscadas 
F(x,y ) = C, expresadas en forma implícita. 

(Lógicamente, para encontrar F podría haberse seguido el proceso anterior cambiando 
el orden en el que se usa P y Q, partiendo de ^ = Q y posteriormente utilizando ^ = P. 
Asimismo, integrando estas dos expresiones e igualándolas, muchas veces basta una simple 
inspección para determinar F .) 


RECETA 4. Ecuaciones exactas. 

Son las de la forma 

P(x, y) dx + Q(x, y) dy — 0, 

es decir, y' = , que cumplen P y — Q x . Se busca una 

función F(x,y) tal que dF — u — Pdx + Qdy, y la solución de la E. D. 
es F(x,y) = C (siendo C constante). 


Ejemplo 4. Resolver 


3 y + e x + (3x + eos y) y' = 0. 


Si ponemos la ecuación en la forma Pdx + Qdy = 0 con P(x,y) = 3y + e x y 
Q(x,y) = 3x + eos y, es claro que P y = Q x = 3, luego la E. D. es exacta. Calculemos 
la función potencial F (que nos dará directamente las soluciones F(x,y) — C ). Como 
Fx = 3y + e x , integrando respecto de x. F(x, y) = 3 yx + e x + <p(y). Derivando respecto de 
y e igualando a Q queda 3x + <p'(y) — 3x + eos y, es decir, ¡p'(y) — eos y, de donde basta 
tomar <p(y) — sen y, y por tanto F(x,y ) = 3 yx + e x + sen y. Así, la solución de la E. D. 
viene dada, implícitamente, por 3yx + e x + sen y — C. (Nótese que al integrar <p'(y) — eos y 
no hace falta poner la constante de integración íp(y) — sen y + C¡ ya que, eu ese caso, una 
de las dos constantes de la solución 3 yx + e x +sen y + C\ — C sería claramente superflua.) 
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Métodos clásicos de resolución de E. D. O. 


APARTADO 4'. 

Reducibles a exactas: Factores integrantes. 

Si tenemos una ecuación P(x, y) dx + Q(x, y) dy — 0 que no es exacta, una idea para 
intentar resolverla sería tratar de encontrar alguna función y(x. y) no idénticamente nula 
tal que 

y(x, y)P{x, y) dx + y(x, y)Q(x , y) dy — 0 

sea exacta. Como esta ecuación es equivalente a la de partida, sus soluciones y las de 
P(x, y) dx + Q(x, y) dy — 0 serán las mismas. 

Desgraciadamente, no hay ningún procedimiento general que permita encontrar 
factores integrantes. Sin embargo, sí que es posible hacerlo, y de manera sencilla, en dos 
casos: 

4'.1. Existencia de factor integrante de la forma y(x). Queremos que / i{x)P{x , y) dx + 
y(x)Q(x, y) dy = 0 sea exacta, esto es, 

d d 

— (n(x)P(x,y)) = — (n(x)Q(x,y)). 
oy ox 

Derivando, y(x)P y (x,y) = y'(x)Q(x, y) + y(x)Q x (x, y), o sea, y(x)(P y (x,y) -Q x (x,y)) = 
y'(x)Q(x,y). Para que esto tenga sentido, 

lAx} _ P y (x,y) - Q x (x,y) 
y(x) Q{x,y ) 

tiene que resultar ser una función que dependa exclusivamente de x, que denotamos h(x). 
Cuando éste es el caso, es claro que la función // que satisface la relación anterior es 

fi(x) — exp ^ j h{x) dx 

con lo cual hemos encontrado el factor integrante buscado. 

4'. 2. Existencia de factor integrante de la forma y(y). Repitiendo el proceso anterior 
buscamos -¡P{y{y)P(x,y)) = j%(l¿(y)Q(x,y)), es decir, y' (y)P + y(y)P y = y(y)Q x , y por 

tanto í j^y — ® x p Py , que tiene que ser función sólo de y, que denotamos h(y). En estas 
condiciones, el factor integrante es y (y) = exp (f h(y) dy). 

4'. 3. Aparte de los casos anteriormente tratados, para algunos tipos de problemas se 
puede intentar encontrar factores integrantes imponiendo a y(x, y) condiciones restrictivas 
de muy diverso tipo. Por ejemplo, exigiendo que sea de la forma y(x, y) = x ot y /3 con a 
y / 3 constantes a determinar, que sea y(x + y), o y(xy), etc. Para estudiar estos casos, 
lo que hay que hacer siempre es igualar yP ) = -§¿{yQ) e intentar resolver la nueva 
ecuación que aparece, teniendo en cuenta, sobre todo, si tiene sentido. Por ser generalmente 
procedimientos bastante particulares, no comentaremos aquí nada más sobre ellos. 
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RECETA 4'. Reducibles a exactas: Factores integrantes. 

Si P(x,y)dx + Q(x,y)dy = 0 no es exacta, podemos intentar 
encontrar p(x, y) tal que 

p(x, y)P(x, y) dx + p(x, y)Q(x, y) dy = 0 

sea exacta. 

4/.1. Existencia de factor integrante de la forma //(.x). Ocurre 
cuando Fy q^ x = h(x), tomándose p(x) — exp (f h(x) dx). 

4'. 2. Existencia de factor integrante de la forma //(y) . Ocurre 
cuando p Píl = h(y ), tomándose p(y) — exp (f h(y) dy). 

4'. 3. Otras expresiones restrictivas para p(x,y). 


Ejemplo 4' . Resolver 

( 2x 2 + y) dx + {x 2 y — x) dy — 0. 

En este caso, P(x,y) = 2x 2 + y y Q(x,y) — x 2 y — x. Esta ecuación no es exacta ya 
que P y — 1 y Q x = 2 xy — 1. Para intentar encontrar un factor integrante se calcula 

P y ~ Qx _ 1 - (2 xy - 1) _ 2(1 - xy) _ -2 
Q x 2 y — x —x{l — xy) x 

Ya que se obtiene una expresión que depende sólo de x. podemos asegurar que existe 
un factor integrante dado por la fórmula p(x) — exp (f ^ dx) — x~ 2 . Entonces, si 
multiplicamos la E. D. por p(x) ~ x~ 2 se obtiene la ecuación exacta (2 + yx~ 2 ) dx + (y — 
x -1 ) dy = 0. Por el método usual, encontramos que la función F tal que F x — 2 + yx~ 2 y 
F y — y — x _1 es F(x,y) = 2x — yx~ l + \y 2 ■ Por tanto, la solución de la E. D. exacta, y 
también de la de partida, resulta ser 2x — yx~ x + \y 2 = C. 


APARTADO 5. 


Ecuaciones lineales de primer orden. 

Dada la ecuación 


y ' + a(x)y = b(x), 


vamos a explicar cómo resolverla por tres métodos distintos: 


(i) Encontrar un factor integrante de la forma //. ( x ) . Para ello, si la ponemos en la 
forma (a(x)y — b(x)) dx + dy — 0 y denotamos P(x,y) — a(x)y — b(x) y Q(x,y) — 1, se 
tiene — = a(x). Por tanto, según liemos visto anteriormente, la E. D. tiene el factor 


integrante p(x) = exp (f a(x)dx). Así, multiplicando por p(x), la ecuación 


exp 


a(x) dx ) (a(x)y — b(x)) dx + exp 


o(x) dx ) dy — 0 
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tiene que ser exacta. Ahora, bastará encontrar la función potencial F con lo que la ecuación 
anterior podrá ponerse dF = 0 y su solución será F(x,y ) = C. Busquemos F: 


Como F y = exp (f a(x)dx), tendremos F — yexp(f a(x) dx) + <p(x). Por otra parte, 
derivando esta F respecto de a; y usando que F x = exp (f a(x ) dx)(a(x)y — b(x)) llegamos 
a 


F x = y exp ( / a(x) dx ) a{x) + <p'(x) = exp ( / a(x) dx ) (a(x)y — b(x )) , 


de donde —b(x)exp(f a(x) dx) — <p'(x). Integrando, <p(x) — — f b(x) exp(f a(x) dx) dx, 
luego F = yexp(f a(x) dx) — f b(x) exp (f a(x) dx) dx y la solución de la ecuación exacta 
(y de la lineal de partida) es, expresada en forma implícita, 


V exp 


a(x) dx 


j b(x ) exp 


a(x) dx ) dx 


C. 


Sin más que despejar y, tenemos que la solución de la ecuación lineal resulta ser 


y = exp 


a(x) dx 


b(x) exp ( / a(x) dx) dx + C 


( ii ) Un segundo método de resolución se basa en resolver previamente la ecuación 
lineal homogénea asociada y' + a(x)y — 0. Esta ecuación es de variables separadas, pues 
puede ponerse ^ = —a(x) dx; su solución es y — C ex p(— f ci(x) dx). 

Apliquemos ahora el método de variación de las constantes , esto es, consideremos 


y — C(x) exp 


a(x) dx 


y vamos a ver cómo debe ser C(x) para que se verifique y' + a(x)y = b(x). Derivando, 
y' = C'(x) exp(— J a(x) dx) — C(x)a(x) exp(— f a(x) dx) y, sustituyendo en la ecuación 
lineal, 


C'(x) exp ( — / a(x) dx ) — C(x)a(x) exp 


a(x) dx 


+ a(x)C(x) exp ( — / a(x) dx ) = b(x). 


Dos de los sumandos anteriores se cancelan, de donde C'(x) = b(x)exp(f a(x) dx) e, 
integrando, 

C(x) = í b(x) exp f í a(x) dx \ + C. 


Así, hemos llegado a la misma expresión para las soluciones que la que encontramos por 
el método anterior. 




E. D. explícitas de primer orden 


19 


(m) El tercer procedimiento de resolución parte de suponer que liemos encontrado, 
de alguna forma, una solución particular y p (x) de la E. D. lineal. Entonces, la solución 
general de la lineal es y p más la solución general de la lineal homogénea asociada, es decir, 


V = V P + C'exp 


a(x ) dx 


es solución para todo C 6 1. La justificación de este hecho es sencilla. En efecto, basta 
comprobar que, si y p es solución de y' + a(x)y = b(x ) y y lo es de y' + a(x)y — 0, entonces 
y + y p es solución de y' + a(x)y — b(x), lo cual es claramente cierto: 

(y + y P Y + a(x)(y + y p ) = (; y ' + a(x)y) + (y p + a(x)y p ) = 0 + b(x) = b(x). 

(iv) El último método de resolución de ecuaciones lineales que describimos consiste 
en efectuar una descomposición y(x) — u(x)v(x) adecuada. Tomando y de esa forma, si 
derivamos, y' = u'v + uv' , con lo cual, al sustituir en la ecuación, u'v + uv' + a(x)uv = b(x). 
Sacando u factor común, podemos escribir la expresión anterior como u'v + (v' + a(x)v)u = 
b(x). Vamos ahora a elegir v de tal forma que se anule el coeficiente de u, es decir, 
que satisfaga v' + a(x)v = 0. Esta es una E. D. en variables separadas; resolviéndola, 
— = —a(x), con lo cual basta tomar logu — — f a(x) dx, es decir, 


v(x) = exp 


— / a(x) dx 


Con v(x) esa función, la ecuación queda ahora u'v = b(x), de donde vi = b(x)v 1 , es decir, 
u'(x) = b(x) exp (f a(x) dx). Integrando, 


u(x) 


b(x) exp 


ci(x) dx 


dx + C. 


Sin más que recomponer y = u{x)v{x), con este procedimiento de nuevo encontramos la 
misma expresión para las soluciones de la E. D. lineal de partida. 

Para concluir, comentemos una vez más que no hace falta recordar la fórmula que 
hemos obtenido para las soluciones de la ecuación lineal, sino que basta seguir en cada 
problema alguno de los procesos descritos. Generalmente, en opinión del que escribe, los 
que suelen conducir a la solución por un procedimiento más corto suelen ser el segundo y 
el cuarto. El tercero tiene, sobre todo, gran importancia teórica. Además, merece la pena 
destacar que el segundo y el tercero tienen su paralelismo a la hora de resolver ecuaciones 
lineales de orden superior (como veremos más adelante), mientras que los otros dos sólo se 
aplican a las de primer orden. 

Ejemplo 5. Resolver 


por los cuatro métodos descritos. 


2 xy' — 3 y — 4x 2 
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RECETA 5. Ecuaciones lineales de primer orden. 

Son de la forma 

y' + a(x)y = b(x). 

Hay tres métodos de resolución: (i) Encontrar un factor integrante 
de la forma y(x). ( ii ) Resolver la ecuación lineal homogénea asociada 
y' + a(x)y — 0 (que es de variables separadas), cuya solución es 
y — (7exp(— f a(x) dx ), y usar el método de variación de las constantes 
(esto es, cambiar C por C(x) en la expresión anterior y sustituir 
en la ecuación lineal), (iii) Encontrar de alguna forma una solución 
particular y p (x), con lo cual la solución general de la lineal es y p 
más la solución general de la homogénea asociada, (iv) Descomponer 
y(x) — u(x)v(x), sustituir en la lineal, e igualar a 0 el coeficiente de 
u, resolviendo la ecuación que aparece (v* + a(x)v = 0, que es de 
variables separadas); tras esto, queda una ecuación en u(x) de variables 
separadas. 

De cualquier modo se obtiene que la solución general de la E. D. 
lineal es 


y = exp 


a(x) dx 


b(x) exp ( / a(x) dx) dx + C 


(i) Tenemos la ecuación lineal y ' — y = 2x, es decir, dy = (y- y + 2 x) dx, que 
es de la forma P dx + Q dy = 0 con P(x,y) = ^ y — 2x y Q(x,y) = 1. Como 

hy q ^ x = existe el factor integrante y{x) — exp (f dx) — x~ 3 ^ 2 . Así, la 

ecuación (P^-x~ 5 / 2 y — 2a: -1 / 2 ) dx + a: -3 / 2 dy = 0 es exacta. La función potencial F debe 
cumplir F y = a; -3 / 2 , luego F — x~ 3 ! 2 y + y>{x). Por otra parte, F x = z ^-x~ 5 / 2 y — 2a: -1 / 2 = 

- 3 - 5/2 


y + de donde ip(x) = —2 f x 1//2 dx — —4a: 1 / 2 . Por tanto, F{x,y) — 

a: -3 / 2 y — 4a; 1 / 2 , y la solución de la E. D. es a; -3 / 2 y — 4a: 1 / 2 — C, o sea y — Cx 3 / 2 + 4a: 2 , 
CgR. 


(ii) La lineal homogénea asociada es y ’ — ^ y = 0. Podemos ponerla como ^ 

de variables separadas, cuya solución es log y = |loga: + K , es decir, y — Cx 3 / 2 . 
Empleemos ahora el método de variación de las constantes, para lo cual tomamos y = 
C(x)x 3 / 2 . Derivando, y ’ — O (x)x 3 / 2 + |C'(a:)a: 1 / 2 y, sustituyendo en E. D. de partida, 
2 x(C’ (x)x 3 ¡ 2 + |C'(a:)a: 1 / 2 ) — 3 C(x)x 3 ^ 2 = 4a: 2 , esto es, C'(x) = 2a: -1 / 2 . Integrando, 
C(x) = 4a: 1 / 2 + Ki luego la solución de la lineal es y = (4a: 1 / 2 + K i)a: 3 / 2 . Si empleamos 
de nuevo C para denotar la constante, y = Cx 3 / 2 + 4a: 2 , la misma expresión que hemos 
encontrado en (i). 

(iii) Primero, tratemos de hallar una solución particular de la ecuación lineal. Lo más 
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sencillo es intentar probar si existe alguna solución polinómica. Esto sólo sería posible 
con un polinomio de segundo grado, pues en otro caso no podrían cancelarse nunca todos 
los sumandos. Por tanto, vamos a tantear con polinomios de la forma y = ax 2 + bx + c. 
Derivando, y' = 2ax + by, sustituyendo en la ecuación, 2x(2ax + b) — 3(ax 2 + bx + c) = 4a; 2 . 
Si igualamos los términos del mismo grado encontramos que esto se cumple con a = 4, 
b — c = 0. Así pues, una solución particular de la lineal es y — 4a; 2 . Por otra parte, 
y tal como ya hemos calculado en (n), la solución general de la homogénea asociada 
es y — Cx 3 / 2 . En definitiva, de nuevo tenemos que la solución general de la lineal es 
y = Cx 3 / 2 + 4a; 2 . 

(iv) Para resolver la ecuación lineal, descomponemos y — u(x)v(x), de donde y' = 
u'v + uv' . Sustituyendo, 2 x(u'v + uv') — 3 uv — 4x 2 , es decir, 2 xu'v + (2xv' — 3v)u = 4x 2 . 
Si igualamos a 0 el coeficiente de u queda 2xv — 3v — 0, o sea integrando, 

log v = |loga;, luego v = a; 3 / 2 . Con esto, la ecuación queda 2 xu'x 3 / 2 — 4a; 2 , es decir, 
v! = 2a; -1 / 2 , cuya solución es u — 4a; 1 / 2 + C. Sin más que recomponer y = uv obtenemos 
y = (4a; 1 / 2 + (7)a; 3 / 2 , la misma solución para la lineal que con los otros tres métodos. 

APARTADO 5'. 

Ecuación de Bernoulli. 

Consideremos 

y' + a(x)y + b(x)y a = 0. 

Es claro que, si a — 0, la E. D. anterior es lineal y, si a — 1, es de variables separadas. En 
otro caso, veamos cómo resolverla: 

Efectuemos el cambio de función y x ~ a — z , para el cual (1 — a)y~ a y' = z' , es decir, 
~^y' — jzt:- Sustituyendo en -^y' + a(x)y 1 ~ a + b{x) = 0 tenemos + a(x)z + b(x) — 0, 
con lo que hemos transformado la ecuación de Bernoulli en la E. D. lineal z'+(l— a)a(x)z — 
(i a — l)b(x). 

Puede seguirse un segundo método de resolución sin más que aplicar un procedimiento 
análogo al último de los que hemos visto para ecuaciones lineales. Partimos de la 
descoposición y(x) = u(x)v(x). Derivando, y ’ = u'v + uv' y, sustituyendo en la E. D., 
u'v + uv' + a(x)uv + b(x)u a v a — 0, que escribimos como u'v + (v' + a(x)v)u + b(x)u a v a — 0. 
Ahora, igualamos a 0 el coeficiente de u. con lo cual tenemos v' + a(x)v — 0, que es una 
ecuación en v(x) de variables separadas; resolviéndola determinamos v(x). Con esta v, 
la ecuación de la que partíamos ha quedado u'v + b(x)u a v a — 0, que es una ecuación 
en u(x) de variables separadas. Resolviéndola encontramos u(x), con lo que ya tenemos 
completamente solucionada la ecuación de Bernoulli. 

Ejemplo 5'. Resolver 

xy' + 2 y + x 5 y 3 e x = 0. 

Puesta en la forma y' + f.y + x A y 3 e x — 0, es claro que la ecuación es de Bernoulli 
con o = 3. Hacemos el cambio de función y~ 2 — z , para el cual z' = —2 y~ 3 y' : es decir, 
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RECETA 5'. Ecuación de Bernoulli. 

Es de la forma 

y' + a(x)y + b(x)y a = 0. 

Si a — 0 es lineal, y si a = 1, de variables separadas. En otro caso, se 
hace el cambio de función y 1-a = z, con lo que la E. D. de Bernoulli 
se transforma en una lineal. Un segundo método de resolución es el 
siguiente: se descompone y(x) — u(x)v(x ) y se sustituye en la E. D., se 
iguala a 0 el coeficiente de u (queda v' + a(x)v = 0, que es de variables 
separadas), lo que nos lleva a determinar v. apareciendo ahora una 
ecuación en u(x) de variables separadas. 


^3 = —y. Si sustituimos esto en p- + fy 2 + x 4 e x — 0 queda z' — — 2x 4 e x , que es 

lineal con a(x) — y b(x) = 2x 4 e x . 

Tal como sabemos, la solución de la ecuación lineal es z — exp (— f a(x ) dx) x 
[/ b(x) exp (f a(x ) dx) dx + C] . Como f a(x) dx — f ^ dx — —4 loga: = loga: -4 , se sigue 
z = x 4 (j 2x x t dx + C) = x 4 (2e x + C ). Por tanto, la solución de la E. D. de Bernoulli es 
y~ 2 — x 4 (2e x + C), es decir, y — x~ 2 (2e x + C) -1 / 2 . 

Resolvamos ahora la ecuación por el segundo procedimiento explicado. Tomamos 
y = uv, luego y' = u'v + uv'\ sustituyendo, x[u'v + uv') + 2 uv + x 5 u 3 v 3 e x — 0, es 
decir, xu'v + (xv' + 2 v)u + x 5 u 3 v 3 e x — 0. Elijamos v tal que xv' + 2v — 0, esto es, 
y = — lo cual se consigue con logn = —2 loga:, o sea, v = x~ 2 . Así, tenemos la nueva 
ecuación xu' + x 5 u 3 x~ 6 e x — 0, que puede escribirse como —(y = e x dx. Su solución es 
2^2 — e x + C / 2, es decir, u = + C) -1 / 2 . Con todo esto, la solución de la ecuación de 

Bernoulli de la que partíamos es y — x~ 2 (2e x + C) -1 / 2 . Lógicamente, hemos obtenido la 
misma que mediante el primer método. 


APARTADO 5". 

Ecuación de Riccati. 

Supongamos que tenemos la E. D. 

y' + a(x)y + b{x)y 2 = c{x). 

Para resolverla, tenemos que haber encontrado previamente una solución particular y p (x). 
Si éste es el caso, efectuamos el cambio de función y = y p + z, con lo cual y' — y' p + z' y, 
sustituyendo, y p +z r +a(x)(y p +z)+b(x)(y 2 +2y p z+z 2 ) = c(x). Como y' p +a(x)y p +b(x)y 2 — 
c(x) por ser y p solución particular, esa expresión queda z' + a(x)z + b(x)z 2 + 2b(x)zy p — 0, 
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es decir, z' + [a(x) + 2 b(x)y p (x)\z + b(x)z 2 — O, que es una E. D. de Bernoulli con a — 2. 
(Nótese ahora que, como el cambio z — u~ l en la de Bernoulli reduciría ésta a una lineal, 
el cambio directo y = y p + - transforma la ecuación de Riccati en una lineal de un solo 
paso.) 

Existen diversos resultados interesantes sobre la ecuación de Riccati. Por ejemplo, si se 
conocen dos soluciones y\ e y 2 suyas, el cambio de función y = H e va a que se pueda 

resolver la ecuación con una sola cuadratura: v = C exp(f b(x)(y 2 (x) — y\(x))dx). Pero 
desafortunadamente, no hay ningún procedimiento general que permita obtener alguna 
solución particular, por lo que muchas veces las E. D. de Riccati 110 resultan integrables 
en términos de funciones elementales. A este respecto, resulta apropiado citar la ecuación 
especial de Riccati y' + by 2 = cx m con b, c G M \ {0}; se sabe que puede resolverse en 
términos finitos si y sólo si m = — 2 o m = 2A . u ''. para algún entero k. 


RECETA 5". Ecuación de Riccati. 

Es de la forma 

y' + a(x)y + b(x)y 2 = c(x). 

El método requiere haber encontrado previamente una solución parti- 
cular y p (x). Si éste es el caso, haciendo el cambio de función y = y p + z, 
la E. D. de Riccati se reduce a una de Bernoulli con a — 2. 


Ejemplo 5". Resolver 


y' + y 2 — x 2 — 2x. 


Un primer vistazo a la ecuación nos induce a pensar que, posiblemente, tenga como 
solución un polinomio de primer grado y = ax + b. Sustituyendo, a + (ax + b) 2 = x 2 — 2x, 
lo que, efectivamente, se satisface para a — — 1 y 6 = 1. Por tanto, una solución particular 
es y p = — x + 1. Si efectuamos el cambio de variable y = —x + 1 + z (luego y' = —1 + z') 
y sustituimos en la E. D., obtenemos, tras simplificar, z' + (2 — 2 x)z + z 2 — 0, ecuación 
de Bernoulli con a = 2. Haciendo el cambio 2 = u~ x (de donde u' — Etz'), al sustituir en 
Erz' + 2x ~ 2 —1 = 0 aparece v! + (2x — 2)it = 1, que es una E. D. lineal con a(x ) — 2x — 2 
y b(x) — 1. La solución de esta última es 


u — exp ( f (2 — 2x) dx) 


exp (/( 2 x — 2) dx) dx + C 


= e 


2x—x z 


^x ¿ —2x 


dx + C 


Deshaciendo los cambios, u — ^ , luego 


y + x — 1 


— e 


2x—x Á 


^x 2 —2x 


dx + C 
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es la solución de la E. D. de Riccati. (De aquí puede despejarse y si se desea, obteniéndose 
la solución explícitamente, pero la integral que aparece no puede resolverse en términos de 
funciones elementales.) 

APARTADO 6. 

Sustituciones. 

Supongamos que tenemos la E. D. 

y' = f(x,y) 

que tiene un aspecto diferente a cualquiera de las que ya se lian estudiado. En ocasiones 
puede ocurrir que, sustituyendo una parte de la ecuación por una nueva variable, la E. D. 
aparentemente difícil de la que partíamos se transforme en una que puede resolverse con 
facilidad. Esto es, en esencia, encontrar un cambio de variables inteligente. Aunque no 
pueden darse regias fijas sobre qué sustituciones usar, si es que hay alguna sustitución 
posible, vale la pena intentar algo cuando no se nos ocurre otro camino. Muchas veces, una 
sustitución muy simple puede ser suficiente. 


RECETA 6. Sustituciones. 

Cuando tenemos una E. D. 

y' = f(x,y) 

que no responde a alguno de los tipos estudiados hasta ahora, a veces 
una sustitución (en esencia, un cambio de variable) más o menos 
ingeniosa transforma la ecuación en una reconocible. Lógicamente, no 
puede darse una regla general pero, en todo caso, merece la pena 
intentar algo. 


Ejemplo 6. Resolver 


dy 

x- y = 

dx 


„.3 

2—e y/x . 

y 


La ecuación no es separable, ni homogénea, ni exacta, ni admite factores integrantes 
fj,(x) o n(y), ni, finalmente, es lineal, de Bernoulli o de Riccati. Sin embargo, si la 
ponemos en la forma xy ~ y — 2-e y l x , nos damos cuenta rápidamente que la derivada 

x y 

del exponente y/x es el primer miembro y 2 y . Así, al efectuar el cambio de función 
dado por la sustitución u — la ecuación queda transformada en v! — 2u~ 1 e u , es 
decir, ue~ u du — 2 dx, ecuación en variables separadas. Integrando por partes resulta 
—ue~ u — e~ u + C — 2x, de donde u + 1 = {C — 2x)e u . Volviendo a las variables iniciales, 
^ + 1 = (C — 2x')e y l x , y por tanto la solución de la E. D. de partida es, expresada en forma 
implícita, y + x — x(C — 2x)e y / x . 




ECUACIONES EN LAS QUE 

LA DERIVADA APARECE IMPLÍCITAMENTE 


Son las que tienen la forma 


F(x,y,y') = 0 . 


APARTADO 7. 

F algebraica en y' de grado n. 

Supongamos que tenemos 

WT + ai(x,y){y') n ~ 1 H f a n - 1 (x,y)y' + a n (x,y ) = 0. 

Si pensamos en la expresión anterior como en un polinomio en y’ de grado n igualado 
a cero, y logramos resolver la expresión algebraica, obtenemos las raíces y ’ = fi(x,y), 
i = 1, 2, . . . , n. Es decir, 

{y' ~ fi(x, y)){y' - f 2 (x, y)) ■ ■ ■ ( y ' - f n (x, y)) = 0. 

Por lo tanto, las soluciones de la E. D. de partida serán las de cada una de las nuevas 
ecuaciones diferenciales y' — fi(x, y) = 0 , i — 1 , 2 ,..., n, que habrá que resolver. De esta 
forma obtenemos n familias uniparamétricas de soluciones. 


RECETA 7. F algebraica en y' de grado n. 

Tenemos 

C y') n + y){y') n ~ 1 H f a n _ i(x, y)y' + ci n {x, y) = 0. 

Resolviéndolo como un polinomio en y ’ de grado n igualado a cero 
obtenemos 

(: y' - /íOó y)){y' ~ Í 2 (x, y))--- {y' - fn(x, y)) = 0. 

Por tanto, las soluciones de la E. D. de partida serán las soluciones de 
cada una de las ecuaciones y 1 — fi(x, y) = 0, i = 1 , 2, . . . , n. 
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Ejemplo 7. Resolver 

y 2 ((y') 2 + 1) = 1- 

2 2" 

Despejando ( y ') 2 queda ( y ') 2 = 1- | , cuyas soluciones algebraicas son y' = ± — — — . 

Vamos a resolver estas dos nuevas ecuaciones diferenciales a la vez (ambas son de variables 

separadas). Si las ponemos como ± J y = dx, integrando, — V 2 = x + C. Si 

v 1— 

elevamos al cuadrado ambos términos, podemos expresar conjuntamente las dos familias 
de soluciones como 1 — y 2 = (cc + C') 2 , es decir, (x + C) 2 +y 2 = 1, que son las circunferencias 
con centro en el eje x y radio 1. 

Por último, es evidente que y = 1 e y — — 1 también son soluciones de la E. D., 
aunque no se encuentran entre las que acabamos de hallar. Si atendemos a la interpretación 
geométrica de las ecuaciones diferenciales, es lógico que estas dos rectas sean soluciones, 
ya que son las envolventes de la familia de circunferencias que satisfacen la E. D. 

Obtención de la envolvente de una familia de curvas. 

Recordemos qué es una envolvente: Dada una familia uniparamétrica de curvas 
Wo} ae A en el plano, la envolvente de la familia es una nueva curva ¡p tal que, en cada 
punto de contacto de p> con alguna de las p> a , la tangente de p> y de <p a es la misma. 

• Supongamos en primer lugar que tenemos una familia de curvas expresadas 
implícitamente como F(x,y,C ) = 0 (es decir, para cada valor del parámetro C aparece 
una curva de la familia). Las envolventes se obtienen eliminando C del sistema 

í F(x, y, C) = 0 
\F c {x,y,C) = 0 ’ 

donde con Fq hemos denotado la derivada parcial de F respecto a C. 

En el caso de las circunferencias del ejemplo anterior tendríamos 

í ( x + C) 2 + y 2 = l 
[2(a; + C) = 0 

Despejando en la segunda ecuación, x = —C y, sustituyendo en la primera, y 2 = 1. Es 
decir, las envolventes son las rectas y = ±1. 

• En segundo lugar, supongamos que tenemos la familia de curvas en paramétricas 

{ x = x(t, C ) 

V = y(t,C) 

Las envolventes se obtienen despejando C = C(t) en el jacobiano igualado a cero 

det( XC »)=(> 

V y t ) 
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y sustituyendo en las ecuaciones paramétricas. 

Realmente, con estos procedimientos no sólo aparecen las envolventes, sino también 
los lugares geométricos de puntos singulares, ya sea de puntos de retroceso, de puntos de 
inflexión o de cúspides. (En los puntos de retroceso no hay vector tangente, y surgen en 
el proceso anterior sea cual sea la parametrización de las curvas; en cambio, los puntos 
de inflexión y las cúspides pueden aparecer o no según sea la parametrización usada. Por 
ejemplo, nunca surgen al usar el parámetro arco. Por esta razón, a los primeros se les 
llama puntos singulares esenciales y, a los otros, puntos singulares evitables .) En todo 
caso, después de aplicar estos métodos, hay que comprobar siempre qué es lo que hemos 
encontrado. 


APARTADO 8. 

Ecuación de la forma y — f(x,y'). 

Como procedimiento general para intentar resolver este tipo de ecuaciones, tomamos 
y' — p y derivamos y — f(x, y') respecto de x, con lo cual tenemos 

P = V' = fx + fy< ^ = fx(x,p ) + fy'(x,p)p'. 

(Quizás resulte más sencillo interpretar esto como que derivamos y = f(x,p) respecto de 
x, obteniéndose directamente p — f x + f P p' ■) Cuando / tiene la forma adecuada, la nueva 
ecuación p = f x (x, p) + f p (x, p)p' que hemos encontrado puede ser de alguno de los tipos ya 
estudiados. Si éste es el caso, la resolvemos, obteniendo su solución x — (f>(p, C ). Entonces, 
la solución de la E. D. de partida será 

\x = <¡>(p, C ) 

\y = f(<f>{p,c),p ) ’ 

expresado como una familia de curvas en paramétricas. Puede resultar extraño pensar que 
el parámetro p vale precisamente -r(; pero esto no importa en absoluto, sino que puede 
considerarse una simple curiosidad. 

Sin embargo, no siempre ocurre, para una función / genérica, que el proceso anterior 
conduzca a una ecuación reconocible. Esto sí que sucede en los tres casos que estudiamos 
a continuación: 


8.1. Ecuación 

y = fW)- 


En este caso, como f(x,y') = f{y'), se tiene f x (x,p ) = 0, luego con el proceso descrito 
habremos obtenido la ecuación p = f p (p)p' o, lo que es lo mismo, p = f'(p)p'. Si p ^ 0, 

esto lo podemos poner como dx — ^ ^ dp , que es de variables separadas. Integrándola, 
x = f ^ ^ dp — <j>(p) + C, y las soluciones de y — f(y') serán las curvas 


x = 4>(p) + c 

y = f(p ) 
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Notar, a la vista de su representación paramétrica, que todas ellas tienen la misma forma: 
se diferencian únicamente en un desplazamiento horizontal. 

A partir de p — 0 obtenemos la solución constante y = /(O) (al sustituir en la ecuación 
queda /(O) = /(O)). Gráficamente, esto es una recta envolvente de las demás soluciones. 

8.2. Ecuación de Lagrange: 

y + x<p(y') + t%') = 0. 

Si ponemos y’ — p y derivamos respecto de x queda 

p + píp) + x<p'{p) ^ + V^'O) ^ = 0 

que, si lo preferimos, podemos escribir como (p + <p(p)) dx + xtp'(p) dp + ip'{p) dp — 0. 

• Supongamos en primer lugar que p + p{p) ^ 0. Entonces, podemos dividir por 
p + p{p), con lo que nos aparece la ecuación lineal 

dx <p'(p) ^'{p) 

— x — U 

dp p + <p(p) p + <p(p) 

en la que x actúa como función y p como variable. Si la resolvemos, obtendremos como 
solución x = ó(p. C). Con esto, habremos encontrado para la E. D. de Lagrange las 
soluciones 

í x = <¡>{p, C ) 

\y = -<j*(p, C)<p(p) - v(p) 

en forma de familia de curvas paramétricas. 

• Supongamos ahora que existe algún A para el cual A + p(A) = 0. Formalmente, 
tomamos y' — Ay, sustituyendo en la E. D., y + xp>( A) + , 0(A) = 0, es decir, y — Xx — ^(X). 
Esta recta es, efectivamente, solución de la ecuación de Lagrange, sin más que comprobar 
que la satisface: Xx — ^(X) +x<p(X) +V’(A) = 0 pues A + <p(A) = 0. Estas rectas decimos que 
son las soluciones singulares de la ecuación. Además, recíprocamente, vamos a ver que si 
y — Xx — ip( X ) es solución de la E. D., entonces se cumple que A + <p(A) = 0. Para ello, no 
hay más que utilizar que y — Xx — L(A) satisface la ecuación, con lo cual deberá cumplirse 
0 = Xx — ipiX) + xip(X) + '¡/K A ) = (A + (p(X))x, luego A + p(X) — 0. 

8.3. Ecuación de Clairaut: 


y ~ xy' + ip(y') = 0. 

Este es un caso particular de ecuación de Lagrange con <p(y') — —y’. Aquí, como 
<p(A) + A = 0 para todo A G M, no existen las soluciones que, en la ecuación de Lagrange, 
encontrábamos por el método general; sólo aparecen rectas. En realidad, tenemos ahora 
como soluciones de la E. D. toda la familia de rectas y — Xx — tí’(A), A G M. 
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Vamos a ver que, además, existe una solución singular, la envolvente de este haz de 
rectas. Para ello, planteamos el sistema 

( y = \ X - 'lp(X) 

| 0 = x — ip'(X) 

en el que la segunda ecuación es la primera derivada respecto al parámetro A. Para 
encontrar explícitamente la envolvente habría que despejar A en una de las dos ecuaciones 
y sustituir en la otra, lo cual no siempre es posible. Pero podemos dejar su ecuación 
paramétrica que, obviamente, será 


í x = /(A) 

\y = Xip'(X) — ip(x) ' 

Aparte de la consideración geométrica de que las rectas no tienen puntos de retroceso, 
cúspides o puntos de inflexión, realmente podemos garantizar que esta curva es una 
envolvente y no un lugar geométrico de puntos singulares. En efecto, en cada punto 
(x(X),y(X)) de la curva, esta interseca precisamente a la recta y = Xx — ip(X) del haz; 
y la pendiente de esa recta en el punto de intersección (y en cualquier otro, de hecho), que 
es A, coincide con la de la curva: 

dy dy/dX ip' (X) + X'ip" (X) — ip 1 (X) 
dx dx/dX ip"( A) 

De paso, esto significa que las soluciones rectas son las tangentes de la solución singular. 

Realmente, podemos hallar las soluciones de la ecuación de Clairaut sin necesidad de 

saber que es un caso particular de la ecuación de Lagrange cuyas soluciones son rectas y su 

envolvente. Para ello, tomamos y' = p en la E. D. y derivamos la expresión y— xp+ip(p) — 0 

respecto de x, obteniendo p — p — xp' + tp'{p)p' — 0, es decir, (—x + tp' {p))p' — 0. De 

,2 

aquí, p' = 0 o ip'(p) = x. En primer lugar, si suponemos -r-f = p' — 0, tendremos 
j¡q: — p = X constante; sustituyendo en la E. D. encontramos la familia de soluciones rectas 
y = Xx — ip( A), A G M. En segundo lugar, si suponemos ip'(p) — x , encontramos la solución 
dada paramétricamente por 

{ x = ip\p) 
y = pip'(p) - ip(p ) 

(donde, anecdóticamente, el parámetro p coincide con p — 4|); ésta es la solución singular, 
la envolvente del haz de rectas. 

Ejemplo 8.1. Resolver 

y = (y') 2 + 2(y') 3 . 


Tomando y' — p podemos poner y = p 2 + 2 p 3 . Derivando respecto de x tenemos 
y' = p — 2 pp' + Qp 2 p r — (2 p + 6 p 2 )p' , es decir, p = (2 + 6 p)pp' . Para p ^ 0 , simplificamos 
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RECETA 8. Ecuación de la forma y — f(x,y'). 

En general, se toma y' = p y se deriva la ecuación y — f(x, y') 
respecto de i. Si / tiene la forma adecuada, a la nueva E. D. se le 
puede aplicar alguno de los métodos ya estudiados, procediéndose así 
a su resolución. 

8.1. Cuando 

y = f(y'), 

la ecuación que se obtiene mediante el proceso anterior es de variables 
separadas. 

8.2. Ecuación de Lagrange: 

y + x<p(y') + ^{y') = 0. 

Se reduce a una ecuación lineal con x como función y p como variable. 
Además, para los A tales que A + <p(A) = 0 se obtienen como soluciones 
las rectas y = Xx — ^(A). 

8.3. Ecuación de Clairaut: 

y - xy' + ip(y') = 0. 

Es un caso particular de ecuación de Lagrange en el que sólo aparecen 
rectas (y su envolvente). 


por p, con lo que 1 = (2 + 6 p)p', es decir, dx — (2 + 6 p) dp. Integrando, obtenemos como 
solución la familia de curvas paramétricas 

í x — 2 p + 3 p 2 + C 

\ y — P 2 + 2 p 3 

Además, a partir de p — 0 se tiene la solución singular y — 0 . 

Para concluir el ejemplo, vamos a ver qué obtenemos si intentamos encontrar la 
envolvente de la familia de curvas paramétricas anteriores. Se tiene 

det (2 + 6p 2 p+V) =P(2+6p)=0 <=> p — 0 o p — y.. 

Con p — 0 obtenemos la recta y = 0 que, efectivamente, es solución de la E. D.: es 
la envolvente de las otras soluciones. Pero con p = aparece la recta y = A- que, 
claramente, no es solución. En realidad, en un punto (x, y) de cualquiera de las curvas, el 
vector tangente es (x,y) = (%, -jjr) = (2 + 6p,2p + 6 p 2 ). Así, cuando p = 0 tenemos el 
vector tangente (i, y) = (2, 0); pero cuando p — ^ el vector tangente es (i, y) — (0, 0), es 
decir, no hay vector tangente. Esto quiere decir que la recta y — es el lugar geométrico 
de los puntos de retroceso de la familia de curvas solución. 
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Ejemplo 8.2. Resolver 

y = x + y'- 3 {y') 2 . 

Nos encontramos ante una ecuación de Lagrange y+x(p(y')+'ip(y / ) = 0 con <p(y') — — 1 
y ^W) — ~y' + 3{y') 2 . Para resolverla, tomamos p = y' y, derivando la E. D. respecto de 

x, tenemos p = 1 + p' — 6 pp' , es decir, (p — 1) dx — (1 — 6 p) dp. Para p ^ 1 podemos poner 
dx — 1 ~^ 6 1 p dp — (—6— jpn) dp, ecuación de variables separadas. (Fijarse que, al describir el 
proceso, hemos dicho que tenía que salir aquí una ecuación lineal; por casualidad, lo que nos 
ha aparecido es una lineal homogénea, que siempre es una E. D. de variables separadas.) 
Integrándola, x — —6 p — 51og(p — 1) + C. Sustituyendo esto y y' = p en la expresión 
y = x + y' — 3(y') 2 queda y — — 6p — 5 log(p — 1) + C +p — 3p 2 — —5p — 3p 2 —5 log(p — 1) + C. 
Así pues, las curvas paramétricas 

f x — —6 p — 5 log (p — 1) + C 
| y — —5 p — 3 p 2 — 5 log (p — 1) + C 

son soluciones de la ecuación de Lagrange de la que partíamos. 

Además de éstas, a partir de p = 1 se obtiene, sin más que sustituir literalmente en 
y — x + p — 3 p 2 , la solución singular dada por la recta y — x — 2. 

Ejemplo 8.3. Resolver 

y = y'x- 2 (y') 2 . 

Estamos ante una ecuación de Clairaut y — xy r + ip(y') — 0 con ip(y') = 2 {y') 2 . Hemos 
demostrado que, sin más que tomar y' = A, las rectas y — \x — 2A 2 son solución de la 
E. D. para cada A G M. Calculemos su envolvente: en el sistema 

í y — Xx — 2A 2 
| 0 = x — 4A 

(recordar que la segunda ecuación se obtiene derivando la primera respecto al parámetro A) 

despejamos A = f en la segunda ecuación y, sustituyendo en la primera, encontramos que 

2 

la envolvente es la parábola y — \- 

Por otra parte, si hubiéramos intentado resolver la E. D. directamente, sin utilizar lo 
que hemos demostrado previamente, tendríamos que tomar y' = p, con lo cual y = px — 2 p 2 

y, derivando respecto de x, p — p'x + p — 4 pp', es decir, (x — 4 p)p' — 0. De aquí, p' — 0 o 

i2 i 

x — 4 p. Si suponemos = p = 0, tendremos ^ = p — A constante, luego, sustituyendo 
en la E. D., hubiéramos encontrado como soluciones las rectas y — Xx — 2 A 2 , A e 1. Por 
último, para hallar la solución singular (envolvente de las rectas), si suponemos x = 4 p, 
tendríamos la solución paramétrica 

{ x — 4 p 

y — px — 2 p 2 = 4 p 2 — 2 p 2 = 2 p 2 

lógicamente, sin más que despejar p en la primera expresión y sustituir en la segunda, esta 
curva es la parábola y = \. 
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APARTADO 9. 


Ecuación de la forma x = f(y,y'). 

Estas ecuaciones son similares a las que aparecen en el apartado 8 pero con el papel 
de x e y intercambiado. En general, para intentar resolverlas, tomamos y r — p y derivamos 
x — f(y,y') (o, si se prefiere interpretarlo así, x — f(y,p )) respecto de y (en lugar de 
respecto a x, como hacíamos en el apartado 8). Esto requiere tener un poco más de cuidado, 
además de usar 4^ = , ), = i. Así, 

ay ay / ax p 7 


i 

p 


dx 

dy 


f v (y,p) + f P (y,p) 


dp 

dy' 


Según como sea /, esta nueva ecuación ^ = f y (y,p ) + f p {y,p ) ^ a la que liemos reducido 
la que teníamos responde a alguno de los tipos previamente estudiados, luego podríamos 
resolverla. En cualquier caso, si logramos hacerlo, obtenemos su solución y = @(p- C). Con 
esto, la solución de la ecuación de partida será la familia de curvas paramétricas 


x = f(HP, C),p) 

y = <t>(p, c) 


Pero no podemos asegurar que, para una función / cualquiera, sepamos resolver la 
ecuación intermedia que nos aparece, con lo que no habría forma de continuar el proceso. 
Se pueden estudiar casos similares a los de la forma y = f(x,y r ) en los que se garantiza 
que el método no quedará interrumpido. Los tres casos que merece la pena distinguir son 
los siguientes: 

9.1. Ecuación x — f(y'). 

9.2. Ecuación x + y<p(y r ) + ip{y') = 0. 

9.3. Ecuación x — jjr + ip(y') = 0. 

Se deja al lector que, como ejercicio, describa el método de resolución de estos tres 
tipos de ecuaciones. Simplemente hay que dedicarse a modificar, con cuidado, el proceso 
utilizado en los correspondientes tipos 8.1, 8.2 y 8.3, recordando que ahora hay que derivar 
respecto de y en lugar de respecto de x. 


RECETA 9. Ecuación de la forma x — f(y, y'). 

En general, se toma y' — p y se deriva la ecuación x — f(y,y') 
respecto de y. Según como sea /, la nueva E. D. que así se obtiene es 
ya conocida, procediéndose a su resolución. Se pueden estudiar casos 
similares a los de la forma y — f(x,y'). 
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Ejemplo 9. Resolver 


x = (y') 3 + y'. 


Si tomamos 4^ = y' 


p y derivamos x — p 3 + p respecto de y tenemos 


1 _ 1 

p dy / dx 


dx 


dp dp 




dy 


dy dy 


que podemos poner como dy = (3 p 3 +p) dp, ecuación en variables separadas. Integrándola, 
y — | P 4 + \p 2 + C. Por lo tanto, las soluciones de la ecuación x = (y') 3 + y' son la familia 
de curvas paramétricas 


x — p + p 


y = I -4 


+ bp 2 + C 


APARTADO 10. 

Ecuación de la forma F(y,y') = 0. 

Para intentar encontrar sus soluciones, consideremos la curva F(a, ¡3 ) = 0. El método 
de resolución requiere que hayamos logrado encontrar previamente una representación 
paramétrica de la curva, esto es, a = ip(t) y ¡3 — ^(t) tal que F(ip(t),ip(t)) = 0. Si así lia 
sido, vamos ahora a efectuar el cambio de función y por t mediante y = <p(t), teniendo en 
cuenta que y' — i/r(t). 

Entonces, derivando y — <p{t) respecto de x tenemos y' = p'(t) ^ que, al ser y' — 
puede escribirse como 

m = £, 

ecuación en variables separadas. 

• Si ip(t) no se anula, podemos poner dx = dt e, integrando, x = / dt + C. 
Como consecuencia, la familia de curvas paramétricas 

í x = /yW cfó + C ' 

1 y = <p{t) 

son soluciones de la E. D. de partida. 

• Si existe algún to para el que ^(¿o) = 0, formalmente efectuamos el siguiente 
razonamiento: 0 = p'{t) así que p'{t) = 0 y consecuentemente <p(t) = cte. = p(to), 
lo que nos conduciría a la solución y — <p(to). Y, en efecto, podemos comprobar con 
rigor que esta recta horizontal es solución de la E. D. sin más que sustituir en ella: 
F{y, y') = F(p(t 0 ), 0) = F(<p(t 0 ), i/>(t o)) = 0. 
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RECETA 10. Ecuación de la forma F(y,y') = 0. 

Consideramos la curva F(a, ¡3 ) = 0. Si encontramos una represen- 
tación paramétrica a = <p(t), (3 — ip(t), F((p(t),ij)(t)) = 0, se hace el 
cambio de función y por t mediante y — <p(t), y ' — ip(t). Así, derivando 
y = <p(t) respecto de x. aparece una ecuación en variables separadas. 


Ejemplo 10. Resolver 

y 2/3 + (y ') 2/3 = l. 

Si tomamos y — eos 3 1 , y' — sen 3 f, es claro que (cos 3 f) 2 / 3 + (sen 3 í) 2 / 3 = 1 . 
Derivando y — eos 3 t tenemos y' = —3 eos 2 1 sen t lo que, usando y' = sen 3 t. resulta ser 
sen 3 t = — 3 eos 2 t sen t . 

ax 

Si suponemos, en primer lugar, sení ^ 0, tenemos dx — ~^° 2 t f dt — — 3cot 2 tdt. 
Integrando x = —3 J cot 2 tdt — —3 f (1 + cot 2 t) dt + 3 J dt = 3 cot t + 3t + C. Por tanto, 
las curvas paramétricas 

( x — 3 cot t + 3t + C 
1 y — eos 3 1 

son soluciones de la E. D. de partida. 

Por último, para los t tales que sen 3 t = 0, lo cual ocurre cuando t = kir, k e Z, 
aparecen como soluciones singulares, tal como justificamos al desarrollar el método, las 
rectas horizontales y — cos 3 (/c7r), es decir, y = 1 e y = —1. 





ECUACIONES DIFERENCIALES EN LAS QUE 
SE PUEDE REDUCIR EL ORDEN 


Supongamos que tenemos la E. D. 

F(x, y, y ', . . . , y^) = 0 con n > 1. 

Al contrario de lo que ocurría con ecuaciones de primer orden, no hay muchos métodos 
para encontrar soluciones de ecuaciones diferenciales de orden mayor que uno. Por ejemplo, 
está demostrado que la E. D. lineal general de orden n > 1 no es soluble por cuadraturas. 
Sí que pueden resolverse las ecuaciones lineales de orden n con coeficientes constantes; 
aunque éste es un tema central en el estudio de la teoría de E. D. O., no nos ocuparemos 
de ello en estas notas. 

Realmente, lo único que vamos a ver aquí es una serie de métodos que permiten reducir 
el orden de una E. D. Aplicándolos (sucesivamente si es necesario), podremos llegar a una 
ecuación de primer orden que, con un poco de suerte, se encontrará entre las que ya sabemos 
resolver. 


APARTADO 11. 

Ecuación de la forma F(x, y ( k \ . . . , y^ n l) — 0 con 1 < k < n. 

Tenemos una ecuación en la que no aparece la variable dependiente y (y además, si 
k > 1, tampoco sus derivadas hasta el orden k — 1). Es evidente que el cambio y ( fc ) — z la 
convierte en 

F(x : z,...,z( n ~ k) ) = 0, 

que es una E. D. de orden n — k. Si logramos resolverla, obtenemos que su solución será 
z — 4>(x, Ci, . . . , C n -k), una familia dependiente de n — k constantes. Entonces, al ser 
y(k) _ para encontrar las soluciones de la ecuación original bastará con integrar k veces 
(f)(x, C\, , C n -k ) respecto de x, es decir, 


y = 



(f)(x, C i, , C n -k ) dx dx . . . dx. 


k veces 


fc veces 


lo cual introducirá k nuevas constantes en la solución general. De este modo, dicha solución 
general tendrá la forma 



y= //•/ <Kx, F \ .... . C n —i¿) ^dX-dx^ .^j—dx, A C n —k~ ^\x A ■ ■ ■ A C n —\x A C n . 

k veces 


k veces 
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RECETA 11. Ecuación de la forma F(x, y^ k \ . . . , y — 0. 

Mediante el cambio y ^ = z se convierte en una ecuación de 
orden n — k. 


Ejemplo 11. Resolver 


y" - xy'" + (y '") 3 = 0 . 


Tomemos y" — z, con lo que la ecuación se transforma en z — xz' + (V) 3 = 0, que es 
de Clairaut. Sin más que sustituir z' — A, las soluciones de ésta son las rectas z — Xx — A 3 , 
Ael. Además, también son solución las envolventes. Para calcularlas, en el sistema 

z = Xx — A 3 
0 = x - 3A 2 


despejamos A = iy/r en I a segunda ecuación, con lo que, sustituyendo en la primera, 
encontramos las dos envolventes z — v/f- Una vez resuelta completamente la ecuación 
de Clairaut, hallemos las soluciones de la original: 

Como y" = z. a partir de z — Xx — A 3 , integrando dos veces, 


V = 


Ar 3 A 3 r 2 

(Ax — A 3 ) dx ) dx — — 1- C\x + C-¿. 

6 2 


Cambiando la notación de las constantes, y — K\x :i — 108 Kfx 2 + K¿x + K.> . 

Por último, de z = obtenemos 

±2 3/2 , A , f ( ±2 x 5 ' 2 \ ±4 xV 2 

y = I l / zt~es x dx) dx — / ( + C\ I dx = + c i x + C 2 . 


3^3 


V3V3 5/2 


15^3 7/2 


Es decir, las dos familias de curvas y = + C\x + C 2 y y = — tüs + C\x + C 2 . 


APARTADO 11'. 

Ecuaciones lineales de orden superior. 

Supongamos que tenemos una ecuación 


a. 


i{x)y (n) + a n _i {x)y { ~ n 1} H h a 1 (x)y r + a 0 (x)y = g(x) 


con n > 1. En primer lugar, conviene dejar constancia de que existe una teoría general 
bien establecida sobre ecuaciones diferenciales lineales de cualquier orden. Pero éste es 
un estudio esencialmente teórico: no es posible describir la solución general de una E. D. 
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lineal de orden 2 o superior por medio de cuadraturas. Sin embargo, sí que se calculan las 
soluciones de forma muy satisfactoria cuando nos encontramos ante ecuaciones lineales con 
coeficientes constantes (es decir, cuando a¿(x) = a¡6l,0<í<fi). Pero todo esto queda 
fuera de los objetivos que estamos persiguiendo aquí. Fundamentalmente, en este apartado 
nos preocuparemos sólo de ver cómo una ecuación lineal de orden n > 1 puede reducirse de 
orden, sucesivamente si es necesario, hasta llegar a una ecuación lineal de primer orden, de 
las que liemos estudiado en el apartado 5. En realidad, este método de reducción de orden 
radica en haber encontrado previamente una solución particular de la lineal homogénea 
asociada 

a n (x)y^ n) + a n _i(x)y (n_1) H f a 1 (x)y' + a 0 (x)y = 0; 

si no hemos sabido hallar tal solución particular, el proceso aquí descrito no tiene para 
nosotros ninguna utilidad práctica. 

Supongamos pues que hemos logrado encontrar una solución particular y n (x) de la 
homogénea asociada (más adelante se pondrá de manifiesto por qué resulta adecuado 
introducir una n en la notación de la solución particular). En este caso, efectuamos el 
cambio de función y por z dado por y(x) = y n (x)z(x). Sin más que derivar sucesivamente, 
obtenemos 


y'{x) = y' n (x)z(x) + y n {x)z'{x) 

y"(x) = y"(x)z(x) + 2 y' n {x)z'{x) + y n (x)z"(x) 


y {n \x ) = yW(x)z(x) 4 + (^jy {k \x)z {n k) (x) 4 + y n (x)z (n) (x). 

Sustituyendo en la lineal, y sacando factor común los z^ k \ queda 

a n ynZ {n) + (na n y' n + a n _iy n )^ (n_1) H h {na n y^~ l) -\ h 2 a 2 y' n + aiy n )z' 

+ («n2/i n) H é ai y' n + a 0 y n )z = g{x). 

El coeficiente de z en la expresión anterior es 0 por ser precisamente y n solución de la 
ecuación lineal homogénea. Así, podemos poner 

b n {x)z (<r ^ + 6 n _i(x)z*' n “ 1 ' ) + • • • + bi(x)z' = g(x). 

Esta es una E. D. de las estudiadas en el apartado anterior. Si en ella hacemos z' = u 
aparece 

6 n (x)'u ( ' n_1 ' ) + 6 n _ i{x)u^ n ~ 2 ^ + • • • + bi(x)u = g(x), 
que es una ecuación lineal de orden n — 1. 

Si sabemos resolver esta nueva ecuación (por ejemplo, si n — 1 = 1 y aplicamos 
cualquiera de los métodos de resolución de E. D. lineales de primer orden ya estudiados), 
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tendremos que su solución general será una función u(x) — (j)(x, C i, , C n - 1 ) dependiente 
de n — 1 constantes. Entonces, la solución de la de partida será 


y(x) = y n (x)z(x) = y n (x) / u(x) dx = y n (x) / <¡>{x, C lt . . . , C n _i) dx 


= Vn{x ) / (¡){x,C 1 ,...,C n - 1 )dx + C n y n (x), 


donde C n surge como constante de integración al calcular una primitiva de ó (que aparece 
representada con la misma notación de integral, esperando que este pequeño y habitual 
abuso de notación no introduzca confusión en el lector). Si no conseguimos resolverla en 
general, se puede intentar reducir el orden de nuevo encontrando ahora u n -i(x) solución 
particular de la homogénea asociada. Si este proceso lo logramos hacer el suficiente número 
de veces, llegaremos siempre a una ecuación lineal de orden 1, que sí que sabemos resolver. 

Finalizaremos comentando cómo varios de los resultados estudiados para E. D. lineales 
de primer orden también pueden aplicarse a lineales de orden superior. En primer lugar, es 
inmediato comprobar (tal como hacíamos en (m) del apartado 5) que la solución general 
de la lineal puede expresarse como una solución particular de la lineal más la solución 
general de la lineal homogénea asociada. (Lógicamente, si la lineal ya es homogénea, puede 
tomarse como solución particular la función nula.) 

Por otra parte, el proceso de reducción de orden que hemos utilizado nos muestra, por 
inducción, que la solución general de la lineal de orden n puede expresarse como 

y(x) = y p (x) + Ciyi(x) H b C n - ij/„_i(a:) + C n y n (x), 

donde y p (x) es una solución particular de la lineal. En efecto, para orden n = 1 lo hemos 
probado al tratar las lineales de primer orden. Y, para efectuar el paso de inducción de 
orden n — 1 a n, basta fijarse en el desarrollo del proceso de reducción seguido. En efecto, 
por hipótesis de inducción, la solución de la E. D. en u(x) de orden n — 1 que nos aparecía 
tendrá la forma 

u(x) = <j)(x,Ci , . . . , Cn— i ) = u p (x ) + CiUi(x) H b C n -iU n -i(x). 

Por tanto, la solución de la ecuación de partida de orden n será 


y{x) = y n (x)z(x) = y n {x) / u(x) dx = y n (x) / cj)(x, C lt . . . , C n - 1 ) dx 


= y n (x) / (u p (x) + CiUi(x) H b C n -iu n -i(x)) dx 


= y n {x) ( / u p (x) dx + C\ / ui(x) dx + ■ ■ ■ + C n -i / u n -i(x) dx + C r 


— Un(x) ( Up(x ) + C\U\{x) + ■ ■ ■ + C n -iU n -i(x) + C n ) 
= y P {x) + Ciyi(x) H b C n -iy n ~i(x) + C n y n (x), 
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tal como queríamos comprobar. 

También existe el método de variación de las constantes para E. D. lineales de orden 
n, y está destinado a resolver la ecuación lineal cuando se conoce la solución general de la 
homogénea asociada, y(x) — Ciyi(x) + C 2 y 2 (x) + • • • + C n y n (x). En estas condiciones, se 
busca la solución de la lineal de la forma 

y(x) = C 1 (x)y 1 (x) + C 2 (x)y 2 {x) H b C n (x)y n (x) 

para ciertas funciones Ck (x) a determinar. Para encontrar estas Ck(x). se deriva 
sucesivamente la expresión anterior y se sustituye en la ecuación lineal, igualando a cero 
los coeficientes de los yk(x). Aunque no nos preocuparemos de ello, puede demostrarse 
rigurosamente que esto siempre conduce a la solución buscada. 


RECETA 11'. Ecuaciones lineales de orden superior. 

Son de la forma 

a n (x)y (n) + a n _i (x)y < ' n ~ 1) H b a 1 (x)y' + a 0 (x)y = g(x). 

Si logramos encontrar alguna solución y n (x) de la lineal homogénea 
asociada, el cambio de función y = y n z hace que la lineal se transforme 
en una del tipo anterior, cuyo orden se puede reducir. Así, aparece una 
nueva ecuación lineal, esta vez de orden n — 1. 


Ejemplo 11'. Resolver 

xy" + ( 7x — 1 )y' — 7 y — x 2 e~ 7x 

sabiendo que la lineal homogénea asociada tiene una solución de la forma e ax . 

En primer lugar, determinemos a para que y 2 = e ax sea solución de xy" + ( 7x — 1 )y' — 
7 y = 0. Derivando, y 2 — ae ax , y" — a 2 e ax ; sustituyendo, ci 2 xe ax + (7x — 1 )ae ax — 7e ax = 0, 
es decir, [x(a 2 + 7ct) — a — 7]e ax = 0, lo cual se consigue con a = —7. 

Apliquemos ahora el método de reducción de orden. Tomamos y = y 2 z, con lo cual, 
derivando, y' = y' 2 z + y 2 z ' , y" = y"z + 2y' 2 z' + y 2 z" . Si sustituimos en la lineal, queda 
(y 2 z + 2y 2 z' + y 2 z")x + (7x — l)(y 2 z + y 2 z') —7y 2 z — x 2 e~ 7x . Sacando factor común z", z' 
y z 1 aparece xy 2 z" + [2xy 2 + (7x — 1 )y 2 \z' + [xy 2 + (7x — 7)y 2 — 7y 2 \z — x 2 e~ 7x . Como el 
coeficiente de z coincide con la lineal homogénea, de la que y 2 es solución, esta ecuación, 
tras sustituir y 2 e y 2 , se transforma en xe~ 7x z" + [—14xe~ 7x + (7x — 1 )e~ 7x ]z' — x 2 e~ 7x . 
Simplificando e^ 7 *, podemos poner xz" + (— 7x — 1 )z' — x 2 . 

(Nótese que suele ser más cómodo utilizar el símbolo y 2 en lugar de su expresión 
concreta ya que, de este modo, no hace falta arrastrar tampoco las expresiones de y 2 ni 
y", sino que basta con dar su valor al final de las simplificaciones. En particular, no se 
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necesita calcular y'¿. Por el contrario, si empleamos directamente las expresiones de y 2 , y 2 
e y 2 desde el principio, esto puede ayudarnos a detectar algún error que hayamos podido 
cometer: es fácil que el error nos llevara a un punto sin salida en el que no se podría 
continuar el proceso al no anularse el coeficiente de z.) 

Tenemos ahora una E. D. de segundo orden en z(x) en la que no aparece explícitamente 
Para reducirla de orden, basta tomar z' — u, con lo cual queda xv! + (— 7x — 1 )u — x 2 , 
esto es, u r + (— 7 — j-)u — x , lineal de primer orden. Resolvámosla aplicando directamente la 
fórmula general que hemos deducido en el apartado 5. Así, sin más que tomar a(x) — —7— -j; 
y b(x) — x, calculando previamente 


a(x) dx'j — exp (/ ( 


— 7 -jdxj = exp(—7x — loga;) — x 1 e 


-l-7x 


tendremos que su solución es 


u{x) — xe 7x j xx 1 e 7x dx + C\ = xe 7x (-^e 7x + C\) — -j-x + C\xe' 


Deshaciendo los cambios, 


y(x) = y 2 (x)z(x) — y 2 (x) / u(x) dx = e rx / (yi + Cixe 7x ) dx 


+ Ci xe 7x dx) — e 


+ lCixe 7x — Ci \e 7x dx 


■ - 7x ( ^x 2 + \C lX e 7x - ±C x e 7x + C 2 ) = f±x 2 e~ 7x + C 2 (±x - ¿) + Re- 


cambiando la notación de C\ por 49C'i . la solución general de la lineal de segundo orden 
queda 

y(x) = Tjx 2 e~ 7x + Ci(7x - 1) + C 2 e~ 7x . 


APARTADO 12. 

Ecuación de la forma F(y. y ', . . . , y ( n )) = 0. 

En esta ecuación no aparece explícitamente la variable independiente x. Vamos a ver 
cómo se reduce de orden si efectuamos el cambio y' — p y la transformamos en una nueva 
en la que aparecerán y, p y las derivadas de p respecto de y. (¡Atención!: respecto de y, 
no respecto de x.) Para ello, vamos a ver el proceso a seguir para escribir las derivadas 
sucesivas de y respecto de x en función de y , p. y las derivadas de p respecto de y. Ya 
tenemos — y' = p. Para las siguientes, si denotamos p' — p" — (yf , . . . , y hacemos 
uso de la regla de la cadena, entonces 

,, dy' dp dp dy , 

V = = tic = dyfa =PPl 
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y 


m 


dy 

dx 




dy 

dx 


[p"p + {p') 2 ] p, 


y así sucesivamente. Nótese que en la expresión de cada y (k l sólo aparecen derivadas de p 
hasta el orden k — 1. 

Con esto, sustituyendo en F los valores que hemos encontrado para y' , y " , . . . , y^ n \ 
la E. D. original se transforma en una de la forma 


G 


í dp d n 1 p\ 
V dy 1 '" 1 dy n ~ l ) 


0, 


cuyo orden es n — 1. 

Si logramos resolver esta nueva ecuación encontraremos que, en general, su solución 
será de la forma p = <j)(y, C \, . . . , (7 n _i), dependiente de n — 1 constantes. Ahora, 
devolviendo a p su valor original, es decir, p = VA obtenemos = cf)(y, Ci, . . . , C n - 1 ), que 
es otra ecuación, esta vez de primer orden, que también hay que resolver (lo cual, además, 
añade una nueva constante). Las soluciones de esta última serán, obviamente, las mismas 
que las de la E. D. de partida. 


RECETA 12. Ecuación de la forma F(y 1 y ', . . . , y^) — 0. 

Hacemos el cambio y' — p y transformamos la E. D. en una nueva 
dependiendo de y, p y las derivadas de p respecto de y. Esta es de 
orden n — 1. 


Ejemplo 12. Resolver 


2 yy" = (y') 2 + i. 


Tomando y' = p, y tal como hemos comprobado al explicar el método, se tiene y" — p'p 
donde p' — Entonces, al sustituir en la E. D., ésta se nos transforma en ‘lyp'p — p 2 + 1, 

es decir, que es de variables separadas. Si integramos, log(p 2 + 1) = log(C'iy), 

de donde p 2 + 1 = C\ y. y por tanto p — F^JC\y — 1. Volvemos ahora a restaurar el valor 
p — dfz, con lo cual nos aparece = ±y/C\y — 1. Esta es de nuevo una ecuación en 
variables separadas, que podemos escribir como [Cpy — 1 )~ 1 ^ 2 dy — ±dx. Integrándola, 
y — 1 = ±C\x + C 2 , con lo que ya hemos encontrado las soluciones de la E. D. de 
partida. 
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APARTADO 12'. 

Si la ecuación F(x, y, y ', ... , y — 0 
es tal que, para a y m fijos, F cumple 

F(\x, X m u 0 , V, . . . , A m - n u n ) = A a F(x, u 0 , u u . . . , u n ), 


vamos a ver que, efectuando un cambio tanto de variable independiente como de 
dependiente, esta E. D. se podrá transformar en una del tipo anterior (es decir, en la 
que no aparecerá la variable independiente). Una ecuación de estas características suele 
decirse que es homogénea generalizada de grado a, en la que cada factor x contribuye con 
grado 1 , cada y con grado m, y' con grado m — 1 , y" con m — 2, etcétera. 

Tomamos dos nuevas variables t y z (t la independiente y 2 la dependiente) que 
relacionamos con las originales x e y mediante x — e t , y — e mí z. Con esto, vamos a ver 
cómo escribir y' = , . . . , y en función de t, z y las derivadas de z respecto de t. 

Si denotamos z' = (y, por supuesto, z ^ ), tenemos 


dy 

dx 


dy/dt 

dx/dt 


me mt z + e mt z' 


= e^ rn - 1)t ( z ' + mz). 


A partir de aquí, 

d 2 y d f dy 


dx 2 dx V dx 


d í dy\ dt 
dt \dx ) dx 


d 


(m-l)t , z ' + mz \\ 

dt \ J 

(m - l)e*- m_1 ) í (z / + mz) + e^ m -^\z" + mz') 

[z" + (2 m — 1 )z' + m(m — 1 )z ) ; 


dx/dt 


— e (m-2 )t( y n 


análogamente 

d 3 y d ( d 2 y\ d í d 2 y\ dt 
dx 3 dx \d 2 x J dt \d 2 x J dx 

= 4- ( e (m “ 2)i (¿" + (2 m - 1 )z' + m(m - 1)*)) — i— 
dt V J dx/dt 

= (■ m — 2)e^ m ~ 2 d (z" + (2 m — 1 )z' + m(m — 1 )z) 

+ e^ m - 2)t ( z '" + (2 m - 1 )z" + m(m - l)z')] - [ - 

J e* 

= e ( rn ~ 3 l t [z!" + (3m — 3)z" + (3m 2 — 6 m + 2)z' + m{m — l)(m — 2)z), 


y así sucesivamente. En general, por inducción, es claro que 


c /( fc )y 

dx^) 


= e im - k)t g k ( 


z, z 


Ák) 


)• 
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De este modo, sustituyendo en la E. D. que estamos intentando resolver, obtenemos 

F(e‘, V + mz ), .... e^~ n >‘g n (z, z',..., a>>)) = 0; 

extrayendo X — e 1 , esto resulta ser 

e at F{l, z, [z! + mz ), . . . , g n (z, z ' , . . . , 2 (n) )) = 0. 

Como e at 110 puede anularse, el otro factor tiene que ser cero, luego liemos transformado 
la ecuación de partida en una de la forma 

G(z,z',...,z^) = 0, 

en la que no aparece explícitamente la variable independiente t. Esta ecuación puede 
reducirse de orden aplicando el proceso descrito en el apartado anterior. E 11 cualquier caso, 
si logramos resolverla, su solución será z = (f>{t, C n ); deshaciendo los cambios x — 

e y — e mt z, tendremos que la solución de la E. D. original será y = x m (j)( loga;, C \, . . . , C n ). 


RECETA 12'. Si la ecuación F(x, y, y 1 , . . . , y^) — 0 

es tal que, para a y m fijos, F cumple 

F{ Xx, X m u 0 , X m ~ 1 ui, . . . , X m ~ n u n ) = X a F{x, u 0 ,u 1 , . . . , u n ), 

haciendo el cambio x = e l , y = e mt z la E. D. se transforma en una 
de la forma G(z, z' , , z^) = 0 , a la que se puede aplicar el método 
anterior. 


Ejemplo 12'. Resolver 


4 x 2 y 3 y" 



Analicemos en primer lugar cómo lograr que esta ecuación sea homogénea generalizada 
en la que cada factor x contribuya con grado 1 , y los factores y , y', y" con grados 
m, m — 1 y m — 2 respectivamente. Para ello, cada monomio tiene que ser del mismo 
grado. En el miembro de la derecha, esto se consigue con 2 = 4m, es decir, rn — 
con lo cual x 2 — y 4 es de grado 2. Entonces, el miembro de la izquierda es de grado 
2 + 3m + (m — 2) = 2 + | + (| — 2) = 2, coincidente con el del otro miembro. Así pues, con 
m = la E. D. es homogénea generalizada (de grado 2, aunque esto no tiene importancia). 


Hagamos ahora los cambios x — e 1 , y — e t ! 2 z, y denotemos z' = 4f, z" — 4r#. De este 
modo, 

. dy dy/dt ^e t / 2 z + e t / 2 z' , /9 , 

« =1 = 1 TS = — ? = 6 (2 + ^ 2) ' 
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Sustituyendo en la E. D., 4 e 2t e 3t / 2 z 3 e st / 2 (z" — jz) = e 2t — e 2t z 4 1 de donde, simplificando, 
4 z 3 z" — 1, ecuación en la que no aparece la variable independiente t. Para resolverla, 
tomamos z' = p, con lo que 


.. dz dp dp dz . 

z = — — = — = ——— p p 

dt dt dz dt 

con la notación p' — Así, tenemos la nueva ecuación Az 3 p'p = 1, con variable 
dependiente p e independiente z, que ya es de primer orden. Esta ecuación, que podemos 
poner como 4 pdp — z~ 3 dz , es de variables separadas. Integrando, 2 p 2 — —^z~ 2 +C i, o sea, 
4p 2 = 2Ci — z~ 2 . Como p — z' — hemos llegado ahora a plantear 4(V) 2 = 2Ci — z ~ 2 , 
es decir, 2 z' — ±\J2C\ — z~ 2 1 que podemos escribir como ±2(2Ci — z ^ 2 ) _1 / 2 dz = dt. De 
nuevo integrando, 


t + C 2 = ± 1 2(2Ci - z~ 2 )~ 1/2 dz — A 1 2z(2C 1 z 2 - 1)“ 1/2 dz 
1 (2Ciz 2 — l) 1 / 2 


'2071 


1/2 


= ±^{2C lZ 2 -l) 1 ' 2 . 


Volviendo a las variables iniciales x — e l e y — é^z, aparece ahora loga; + Ci — 
±¿(2 C iy 2 x~ x - l) 1 / 2 , con lo cual, sin más que cambiar la notación de las constantes, 

K i log x + K '2 = ±\j2Kiy 2 x- { — 1. Así pues, elevando al cuadrado y despejando y 2 , las 
soluciones de la E. D. de partida vienen dadas por 

2 (Ki loga; + K 2 ) 2 + 1 
V ~ 2K lX 

donde K \ y K 2 son constantes arbitrarias. 


APARTADO 13. 

Si la ecuación F(x, y, y ', ... , y = 0 
es tal que, para a fijo, F cumple 

F(x, Xu 0 , Xui, . . . , A u n ) = X a F(x, u 0 ,ui,..., u n ), 

vamos a comprobar cómo, a través de un cambio de variable dependiente, el orden se 
puede reducir en uno. La propiedad que caracteriza a la función F suele expresarse 
diciendo que F es homogénea de grado a respecto a los argumentos uq, u±, . . . , u n . Con este 
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lenguaje, la ecuación que estamos intentando resolver es homogénea de grado a respecto 
de y, y ', . . . ,y (n) . 

En primer lugar, es claro que, cuando a > 0, la función constante y — 0 (para 
la cual y' — y" = ■■■ = 0) es solución ya que, efectivamente, F{x, y, y ', . . . , y^) = 
F{x, 0, 0, . . . , 0) = 0 a F(x, 1, 1, . . . , 1) = 0. 

Para hallar las demás soluciones, tomemos una nueva función z dada por y' = yz. 
Obviamente, esto es equivalente a decir y — exp (f zdx), puesto que 



dy 

y 


— zdx 


log y — zdx 



En cualquier caso, si derivamos sucesivamente la expresión y' = yz , obtenemos y" = 
y'z + yz' — yz 2 + yz' = y(z 2 + z'), y’" — y" z + 2 y'z 1 + yz" — y(z 2 + z')z + 2 yzz' + yz" — 
y{z 2 + 3 zz' + z"). etcétera. Es importante destacar el hecho de que siempre aparece una 
relación del tipo y ^ = ygk(z : z ' , . . . , Entonces, sin más que sustituir, la E. D. 

queda 

F{x, y, yz, y(z 2 + z'),..., yg n (z, z' z (n_1) )) = 0, 
de donde, extrayendo A — y, 

y a F(x, 1, z,z 2 + z',..., g n (z, z (n_1) )) = 0, 


y por tanto, como estamos suponiendo que y no es la función nula, el segundo factor habrá 
de ser cero. Claramente, esto puede ponerse en la forma 

G(x, z, z ' , . . . , = 0, 


que es una ecuación de orden n — 1. Si logramos resolverla, tendremos que sus soluciones 
serán z = (f>(x, C \, . . . , C n - 1 ), una familia dependiente de n — 1 constantes. Entonces, las 
soluciones de la E. D. original serán 


y = exp ^ J <j)(x, C i, . . . , CVi-i) dxj , 
lo que, al integrar, introduce la n-ésima constante. Así, puede ponerse 


y = C n exp 


<t>(x, Ci, . . . , C n -i) dx 


RECETA 13. Si la ecuación F(x, y, y', . . . , y^) — 0 
es tal que, para a fijo, F cumple 

F(x, Xu 0 , Xui, . . . , Xu n ) = X a F(x, u Q , u u . . . , u n ), 

entonces el cambio de función dado por y' = yz (es decir, y — 
exp (J zdx)) hace que el orden se reduzca en uno. 
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Ejemplo 13. Resolver 

3 X \{y'f - yy") = y\ 

Esta ecuación es homogénea de grado 2 respecto de y, y ' , y" pues, al sustituir y , y', y" 
por A y, X y', X y", la ecuación queda multiplicada por A 2 . Para resolverla, hacemos el cambio 
y por z dado por y = exp (f z(x) dx ), con lo cual 


y' — z exp ^ J zdx), y" = (z' + z 2 ) exp (^J zdx 


Sustituyendo en la E. D. queda 


3ah 


z 1 exp 2 / zdx) — (V + z 2 )exp 2 zdx 


— exp 2 zdx) , 


de donde, simplificando exp(2 J z dx), se sigue —3 x 2 z' = 1, ecuación en variables separadas. 
Tras ponerla como dz = — |x -2 dx la integramos, con lo que z = \x~ x + C \ . Por último, 


y = exp 


= exp 


z(x) dx) = exp (|x 1 + C\) dx) 
Q log x + C\x + C‘¿) = x 1 / 3 exp((7ix + C 2 ), 


luego la solución de la ecuación original es, sin más que cambiar la notación de las 
constantes, la familia de funciones y = K\x x ^e K2X . 
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APÉNDICE: Métodos básicos 
para calcular integrales indefinidas 


l.° Reglas principales de integración: 

(1) Si F'(x) = f(x), entonces 

I f(%) dx = F(x) 4- C, 
con C una constante arbitraria. 

(2) J af(x) dx — a j f(x) dx. donde a es una constante. 

(3) J \f(x) ± g{x)) dx = J f(x)dx± J g(x)dx. 

(4) Si j f(x) dx — F(x) + C y u — <p(x), entonces 

í f(u ) du — F(u) + (7, 


siendo du = <p'{x) dx. En particular, si F'{x) = f(x), se tiene 

í f{ax + b) dx — —F(ax + b) + C, a ^ 0. 

/ a 


2.° Tabla de integrales inmediatas: 

/ fc+i 

x k dx — \- C, k —1. 

k + 1 

( b ) f — = log |x| + C. 

J x 

(c) í e x dx — e x + C. 
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(d) / a x dx — b C, a> 0. 

3 log a 

(e) í sen x dx — — cosa; + C. 


(f) / cosa ; dx — senx + C. 


eos 2 x 


— / (1 + tg 2 x ) dx — tg x + C. 


LbtAs I / - 2 \ 7 

— = / ( 1 + cotg x) dx = — cotg x + G . 

sen- x J 

dx _ . | ( X 71 \ 

= log sec x + tg x + C = log tg - + - + C. 

eos x V 2 4 / 

dx | „ x 

= log cosec x — cotg a; + G = log tg — + G . 

senx 2 


(k) / sh x dx = ch x + C. 


(1) / chx dx = slix + C. 


— =th x + C. 


dx , „ 

— 2 — = — coth x + G . 

sh x 


dx 1 x 1 x 

5 5- = — arctg — b G = arccotg — b Gi, a 0. 

j + a 2 a a a a 


dx 1 . x — a ^ 

— — 2 = h - lo S — “ + C i a ^ °- 
x 2 — a 2 2a x + a 


dx 1 a + x „ 1 , x ^ 

5 5- = — log + C = - arg th - + G, a / 0. 

- — x 2 2a a — x a a 


/x 2 + a 


/ 3. X 

= log x + \J x 2 + a + C = — arg sh —= + C, a > 0. 

a -v / a 


= = log x + a/x 2 — a + C — — arg cli — — + C, a > 0. 

'x 2 — a a \/a 


/ d,T X X 

= = are sen — b C — — are eos — b Ci, a > 0. 
Va 2 - x 2 a a 
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3. ° Método de sustitución: 

Si en la integral J f(x)dx se efectúa la sustitución (o cambio de variable ) x — ip{t) 
se tiene 

J f{x) dx — J f(<p(t))<p'{t) dt. 

La función cp se procura elegir de tal manera que el segundo miembro de la fórmula anterior 
tenga una forma más adecuada para la integración. 

4. ° Integración por partes: 

Si tenemos dos funciones u = <p(x) y v — ip(x), se verifica 


/ 


udv = uv — 


I 


v du , 


donde du — <p'(x)dx y dv — ^'(x)dx. En este método, u y v hay que buscarlas de tal 
forma que la nueva integral que aparece sea más asequible que la de partida. 




